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1 .  INTRODUCTION 


This  report  summarizes  research  work  performed  at  the 
Westinghouse  SAD  Center  under  ONR  Contract  No.  N(XX)14-85-C-0035 ,  "Atomic 


Line  Spectral  Filters  for  Neodymium  Laser  Receivers,"  for  the  period 


/ 


between  November  26,  1984  and  November  25,  1985.  The  major  efforts 
were:  (1)  to  reassess  the  technical  merits  of  a  TI  atomic  line  filter 
(ALF)  employing  a  low  vapor  density  T1C1  and  (2)  to  measure  the  excited 
state  cesium  trapping  lifetimes  and  to  derive  the  self ^quenching  cross 
sections  of  excited  state  cesium  atoms. 

These  two  separate  efforts  address  alternative  approaches  to  the 


implementation  of  an  ALF  based  on  metastable  T1  atoms,  Tl^  whose 


absorption  wavelength  (535  nm)  matches  the  frequency  doubled  Nd:Bel 
laser.  The  first  approach  uses  photodissociation  to  produce  Tl^ 
directly  from  T1C1.  The  second  approach  involves  resonantly  pumping 
the  Cs  in  a  Tl-Cs  mixture  and  collisionally  transferring  energy  to  the 


T1  to  sustain  the  required  density  of  T1 


5^ 


I, 


<r 


Previous  experimental  work  by  Liu  (1985)  on  the  first  approach 


confirmed  the  general  viability  of  the  scheme,  but  at  the  T1C1  densities 
43 


15  3 

-  6  X  10  /cm  )  the  metastable  lifetime  was  found  to  be 


used  (2  x  10 

controlled  by  collisional  quenching  by  the  T1C1  molecules.  An  important 
conclusion  from  the  work  was  that  a  T1C1  based  filter  would  need  to  be 
operated  at  lower  T1C1  densities.  Consequently,  the  medium  would  be 
optically  thin  to  the  pumping  light,  and  one  might  expect  an  optimum 
density  to  exist.  The  objective  of  the  modelling  calculations  reported 
in  Section  2  was  to  explore  this  concept  and  establish  where  in 
parameter  space  the  optimum  operating  point  is  likely  to  lie.  In  the 
analysis  we  have  attempted  to  include  all  the  processes  which  might  be 
conceivably  important  in  determining  the  sustainable  Tl#  density,  and  to 
track  the  densities  of  all  the  products.  Of  particular  interest  is  the 


"V 


m 
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density  of  ground  state  thalliua,  T1qI  resulting  directly  froa 
photodissociation  or  froa  quenching  of  the  Tl^.  This  species  is 
expected  to  react  with  the  available  Cl  or  CI2  to  ref ora  T1C1.  Under 
many  of  the  conditions  considered  the  primary  reformation  reaction  is 
with  Cl,  at  the  wall.  This  implies  a  degree  of  coverage  of  the  surface 
with  condensed  Tl,  which  must  be  kept  sufficiently  low  that  it  does  not 
affect  the  optical  properties  of  the  windows.  The  coverage  is 
determined  by  the  balance  between  the  rate  of  condensation  of  the  Tl  and 
its  chemical  ablation  by  the  Cl  atoas.  It  will  therefore  be  a 
decreasing  function  of  teaperature,  and  an  increasing  function  of  [T1q] , 
the  steady  state  volume  concentration  of  ground  state  thallium  atoms. 

The  prescribed  density  of  T1C1  sets  a  lower  bound  to  the  temperature, 
and  this  will  probably  be  sufficient  to  maintain  a  suitably  low  coverage 
of  Tl.  An  upper  bound  on  the  required  teaperature  will  be  the  dew  point 
of  Tl  corresponding  to  the  calculated  value  of  [T1q] .  In  any  case  the 
potential  problem  of  Tl  condensation  is  best  mitigated  by  choosing 
conditions  which  do  not  iaply  excessively  large  values  of  [T1q] .  To 
this  end,  we  explore  the  possible  advantage  of  adding  chlorine  to  the 
material  mix,  as  a  method  of  accelerating  the  reformation  of  T1C1  from 
T1q  and  CI2  via  a  gas  phase  reaction. 

The  results  indicate  that  for  T1C1  vapor  densities  below 
13  -3  10  -3 

10  cm  a  Tl  density  of  10  ca  can  be  sustained  by  a  pump  power 

*  2 

density  of  100  mV/cm  .  This  power  could  be  easily  provided  by  a  mercury 
254  nm  resonance  lamp.  The  required  input  data  to  the  model  are, 
however,  incomplete  and  we  recognise  that  the  conclusions  are  therefore 
correspondingly  uncertain.  The  effects  of  order-of -magnitude  changes  in 
the  various  input  parameters  is  explored  in  order  to  convincingly 
identify  the  more  important  processes.  Thus  the  calculations  identify 
the  approximate  region  where  operation  of  such  a  filter  will  be 
optimised,  and  thereby  will  improve  the  focus  of  future  experiments 
bearing  on  this  approach. 

In  Section  3  we  describe  the  design,  performance  and  analysis  of 
an  experiment  to  measure  the  effective  lifetime  of  Cs  atoms  resonantly 


pumped  to  the  ^P^/2  or  ^3/2  8^a^es  pure  cesium.  This  lifetime 
is  important  in  determining  the  pump  power  requirements  of  a  Tl-Cs  ALF, 
and  therefore  the  experiment  described  is  an  essential  first  step  in 
demonstrating  the  feasibility  of  such  a  device.  It  was  anticipated  from 
the  outset  that  the  two  Cs  excited  states  would  be  collisionally  mixed, 
and  that  their  effective  lifetime  would  be  determined  by  the  effects  of 
radiation  trapping.  In  general  the  latter  effect  is  enhanced  by 
increased  density,  up  to  the  region  where  collisional  broadening  of  the 
lines  becomes  dominant.  In  order  to  quantify  this  behavior  the 

11  -3 

experiments  were  performed  over  a  wide  range  of  Cs  densities,  ~10  cm 
14  -3 

to  "2  X  10  cm  In  the  analysis  of  the  data  we  find  it  necessary  to 

include  the  effects  of  the  hyperfine  structure  of  the  lines  on  the 

radiation  trapping,  and  to  invoke  a  surprisingly  large  cross  section  for 

collisional  self  quenching.  The  presence  of  this  unanticipated  loss 

process  causes  the  effective  lifetime  of  the  Cs*  to  have  a  maximum  value 

13  -3 

of  M35  fiB  occurring  at  a  Cs  density  of  “2  X  10  cm  in  the  present 
apparatus.  In  the  absence  of  collisional  quenching  the  lifetime  would 
reach  80  /is.  The  lower  value  achieved  is  nevertheless  adequate  for  a 
viable  Cs-Tl  ALF. 

Section  4  of  this  report  discusses  the  conclusions  to  be  drawn 
from  this  work,  and  the  resulting  recommendations  for  subsequent  work 
directed  towards  the  ultimate  goal  of  demonstrating  a  viable  ALF  based 
on  Tla,  and  providing  the  data  base  required  to  properly  optimise  its 
final  design. 


2.  THE  T1C1  SYSTEM 


A  simplified  energy  level  diagram  of  thallium  chloride  and 

thallium  is  shown  in  Figure  2.1.  An  active  thallium  ALF  with 

dissociative  excitation  of  T1C1  requires  two  photon  absorption 

processes.  Initially,  the  metastable  state  of  T1  is  populated  by  a 

250  nm  pump  source.  The  absorption  of  a  535  nm  photon  from  a  laser 

transmitter  is  ultimately  detected  via  the  emission  of  a  378  nm  photon 
2 

from  the  Tl(7  state.  The  results  of  previous  studies  (Liu  et  al., 

1085)  on  T1  atomic  line  filters  using  T1C1  indicated  that  the  pump  power 
at  250  nm  required  for  photodissociative  excitation  increases  with  the 
T1C1  vapor  density  due  to  a  large  collisional  quenching  cross  section  of 
TlC^P^g)  atoms  by  T1C1  molecules.  In  order  to  reduce  the  pump  power, 
we  have  explored  the  possibility  of  using  very  low  T1C1  vapor  density. 
One  conclusion  from  the  study  was  that  the  optimal  conditions  for  a  T1C1 
based  filter  would  involve  lower  T1C1  vapor  density  from  10  /cm  to 
below  10  /cm  .  Loss  of  metastables  by  diffusion  to  the  walls  could  be 
prevented  by  addition  of  a  buffer  gas.  Starting  from  this  modified 
concept  of  how  to  implement  this  particular  atomic  line  filter  we 
performed  the  modeling  calculations.  We  have  attempted  to  include  in 
the  calculations  all  conceivably  important  processes.  The  required  data 
is,  however,  incomplete  and  we  recognise  from  the  outset  that  our 
conclusions  carry  a  corresponding  uncertainty.  Nevertheless,  it 
fulfills  the  primary  objective  of  this  study,  namely  to  identify  the 
approximate  region  of  parameter  space  where  operation  of  such  a  filter 
will  be  optimized.  Exploration  of  this  region  in  future  experiments 
should  then  be  able  to  remove  the  uncertainties  inherent  in  the  present 
calculations . 


2.1  STBADT  STATB  ANALYSIS 


In  designing  a  thallium  atomic  line  filter  the  quantities  of 
primary  importance  are:  (1)  the  achievable  density  of  the  thallium 


metastables  Tln(=[Tlm] ) ,  (2)  the  required  pumping  photon  flux  density 


(=Np) ,  and  (3)  the  operating  temperature  required  to  keep  all  the  active 


species  in  the  gaseous  state.  In  this  analysis,  all  conceivable 
reactions  which  might  affect  these  quantities  are  taken  into  account. 
They  are  listed  in  Table  2.1  as  Eqs.  (2.1)  to  (2.11)  and  in  each  case 
the  relevant  reaction  cross  section  or  rate  coefficient  is  indicated. 


In  addition,  we  include  the  effects  of  diffusion  of  T1  ,  T1  ,  and  Cl  to 
’  m’  o' 


the  walls  of  the  container,  where  T1  is  assumed  to  be  quenched  with 

m 

probability  p^,  and  T1q  and  Cl  are  assumed  to  stick  with  probabilities 
Pj  and  pg.  Under  steady  state  conditions  T1q  and  Cl  stick  at  equal 
rates,  and  leave  as  T1C1  at  the  same  rate. 

The  complete  set  of  reactions  is  represented  schematically  in 
Figure  2.2.  For  the  purpose  of  this  analysis,  it  is  assumed  that  the 
steady  state  spatial  distributions  of  the  relevant  species  are 
adequately  represented  by  the  fundamental  diffusion  mode  of  the 
container,  modified  by  the  particle  reflection  at  the  physical 
boundaries.  With  this  assumption  we  may,  in  general,  represent  the 


local  diffusion  loss  rate  of  species  n^  by 
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Table  2. 


1  Reactions  involved  in  sustaining  a  T1  density  by 
photo-dissociation  of  T1C1.  See  also  figure  2.2 
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In  these  equations  is  the  diffusion  coefficient  of  particles  n^, 
assumed  to  be  independent  of  position,  and  A^  is  the  effective 
fundamental  mode  diffusion  length  of  the  container.  In  calculating 
account  must  be  taken  of  the  non-sero  particle  density  at  the  walls.  At 
high  densities  of  scatterers,  N,  and  provided  the  sticking  probability 
p^  is  close  to  unity,  then  A^  differs  little  from  the  conventional  value 
Aq.  In  the  present  problem  neither  assumption  is  necessarily  valid.  In 
this  situation  it  is  appropriate  to  employ  a  boundary  condition 
involving  the  concept  of  a  * linear  extrapolation  distance,*  X^,  whose 
inverse  is  the  logarithmic  gradient  of  the  particle  density  at  the  wall: 


1 


n. 

1 


(2.14) 


The  definition  of  X^  in  terms  of  the  diffusion  mean  free  path 
(or  equivalently  the  diffusion  coefficient  D.) ,  and  the  sticking 
probability  p^  has  been  discussed  in  detail  elsewhere  (Chantry,  1985). 
For  the  present  purposes  we  employ  the  empirical  analytical 
approximation 
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+  1  X. 
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(2.15) 


where  Aq  is  the  container  diffusion  length  computed  in  the  conventional 
way,  i.e.  with  the  assumption  that  n^(wall)  =  0,  and  1Q  is  the  container 
volume  to  surface  area  ratio.  For  the  cylindrical  containers  of  present 
interest  these  quantities  are  given  by 


(2.16) 
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(2.18) 


Use  of  Eqs.  (2.15)  through  (2.18)  is  expected  to  giro  values  of 
A.  accurate  to  10%  or  better.  In  view  of  the  other  uncertainties  in  the 

l 

present  analysis  this  degree  of  accuracy  is  considered  acceptable.  Use 
of  Eq.  (2.15)  avoids  the  need  to  solve  iteratively  the  transcendental 
equations  involved  in  a  case  by  case  application  of  Eq.  (2.14). 

In  the  following  analysis  we  shall  use  the  subscripts  i  =  a,  1, 
and  2  to  denote  quantities  representative  of  species  Tl^,  T1q,  and  Cl 
respectively.  For  the  tiae  being  we  represent  their  respective  local 
diffusion  loss  frequencies  by  d^,  d^,  and  d 2>  with  the  understanding 
that  their  values  are  quantitatively  determined  by  Eq.  (2.13),  and 
(2.15)  through  (2.18). 

In  developing  steady  state  solutions  to  Eqs.  (2.1)  to  (2.10)  we 


adopt  the  viewpoint  that  the  pumping  photon  flux  N^,  and  densities 


[T1C1] ,  [Clj] ,  and  [M]  (where  X  represents  an  added  inert  buffer  gas 
molecule)  can  all  be  prescribed.  This  assumption  with  regard  to  Clg 
implies  the  possibility  of  adding  additional  chlorine,  in  excess  of  T1C1 
stoichiometry.  Hence,  equating  production  and  loss  rates  for  T1C1 
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Balancing  the  rates  for  Clj: 


gives 
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Balancing  the  rates  for  T1 
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For  completeness  we  Bay  also  solve  for  dj,  the  diffusion  loss  frequency 
of  Cl,  by  equating  production  and  loss  of  Cl: 
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(2.22) 


We  can  now  proceed  to  calculate  all  the  steady  state  densities 
as  follows: 

1.  Prescribe  molecular  densities  [T1C1] , [Clg]  and  [M] ,  all  the 

relevant  reaction  rate  parameters,  p  and  p« ,  and  the 

m  x 

container  dimensions. 

2.  Compute  dj  using  Eqs.  (2.13),  (2.15)  through  (2.18)  and 
N  =  [M  +  T1C1  +  Cl] . 

3.  Compute  [T1q]  from  Eq.  (2.10). 

4.  Compute  [Cl]  from  Eq.  (2.20). 

5.  Compute  d^  using  Eq.  (2.13),  Eq.  (2.15)  through  (2.18)  and 
N  =  [M  +  T1C1  +  Cl] . 

6.  Compute  [Tl^]  from  Eq.  (2.21). 

7.  Repeat  steps  2  to  6  to  solve  iteratively  for  [Cl],  [T1q] , 
and  [T1J. 

8.  Compute  Pg  by 

a.  computing  d„  from  Eq.  (2.22),  then 

b.  computing  I/A2  from  Eq.  (2.13),  then 

c.  computing  X2  from  Eq.  (2.15),  then 

d.  computing  P2  from  Eq.  (2.18). 

Step  7  is  an  iteration  loop  which  allows  Cl  and  T1q  to  be 

included  as  constituents  of  the  total  neutral  density  N,  which  in 

general  is  taken  to  be  [M]  +  [T1C1]  +  [T1  ]  +  [Cl„]  +  [Cl] .  For  the 

first  iteration  [Cl]  +  [T1  ]  +  [Cl„]  is  set  to  sero. 
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2.2  INPUT  DATA 


2.2.1  Absorption  Cross  Sections 

The  total  absorption  cross  for  T1C1  has  been  measured  as  a 

function  of  wavelength  by  Davidovits  and  Bellisio  (1969) .  Ve  have 

assumed  that  the  process  peaking  at  a  wavelength  of  251  nm  with  a  value 
-17  2 

of  5  x  10  cm  is  entirely  dissociative.  The  branching  of  this 
process  has  been  measured  by  Van  Veen  et  al.  (1981)  at  248  nm,  giving 
71%  into  the  metastable  channel.  Based  on  these  numbers  we  assign  the 
following  values: 

al  =  3.5  x  10‘17  cm2 
*2  =  1.5  X  10"17  cm2 

We  note  however  that,  in  a  paper  published  after  completion  of 
the  present  calculations,  Scott  and  Piper  (1986)  have  reported  evidence 
that  only  11%  of  the  absorption  at  249  nm  leads  to  atomic  fragments. 

They  draw  similar  conclusions  regarding  TlBr  and  Til,  some  of  which 
appear  to  be  irreconcilable  with  the  fluorescence  quantum  efficiency 
measurements  of  Maya  (1978).  If,  according  to  Scott  and  Piper  (1986), 
the  values  of  a ^  and  a g  are  an  order  of  magnitude  smaller  than  indicated 
above,  the  present  conclusions  would  need  to  be  modified,  primarily 
regarding  the  required  pump  power.  The  differences  are  apparent  from 
the  sensitivity  study  described  in  Section  2.3.5  where,  among  other 
things,  the  effect  of  decreasing  these  cross  sections  by  an  order  of 
magnitude  is  explored.  Since  we  are  concerned  here  only  with  optically 
thin  situations  the  total  absorption  cross  section  is  not  of  direct 
significance. 


2.2.2  Quenching  Cross  Sections 

Quenching  of  Tl^ g)  by  a  number  of  gases  has  been  measured 

by  Bellisio  and  Davidovits  (1970) .  Their  value  for  T1C1  is 
-15  2 

3.9  X  10  cm' .  Liu  (1985)  obtained  results  consistent  with  this 

value,  or  possibly  larger.  In  order  to  treat  this  important  reaction 

—14  2 

conservatively  we  assume  a  value  of  <7g  =  1  X  10~  cm  for  most  of  our 
calculations,  and  in  addition  identify  optimised  conditions  when 
ff3  =  4  x  10~15  cm2. 

We  are  not  aware  of  any  published  data  for  reaction  (2.4). 

-14  2 

Gedeon  et  al.  (1971)  have  measured  a  value  of  1.59  X  10~  cm  for  the 

analogous  reaction  Tl^O^g/g)  ♦  ^  *  T1Q  +  ^2‘  therefore  assume  a 

value  of  a4  =  1.5  x  10  14  cm2. 

The  self -quenching  cross  section,  reaction  (2.7),  has  a 

surprisingly  large  cross  section  considering  the  large  energy  defect 

-16  2 

(0.97  eV) .  A  value  of  5.4  X  10  cm  published  by  Pickett  and  Anderson 
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(1969)  is  supported  by  the  value  of  4.5  X  10~  cm  published  by 
Aleksandrov  et  al .  (1978) .  We  therefore  use  a  value  of 
a7  =  5  x  10-16  cm2. 

Aleksandrov  et  al.  (1978)  have  also  measured  quenching  cross 

sections  for  a  range  of  buffer  gases.  Their  values  for  N„  and  Ar  are 

-20  2  -24  2  * 

5  X  10  cm  and  5  x  10  cm  respectively.  Either  gas  would  probably 

be  satisfactory  as  the  buffer  M  in  the  present  application.  We  assume  a 

-19  2 

conservative  value  of  ffg  =  1  x  10  cm  . 

In  the  absence  of  any  data  on  reaction  (2.9)  we  have  assumed  a 
-17  2 

value  of  <7g  =  1  X  10  cm  which  is  probably  conservative  based  on 
values  for  other  atomic  quenchers. 

2 . 2 . S  tef ormati on  Keactions 

We  are  not  aware  of  data  for  reaction  (2.5).  The  analogous 

reaction  T1  +  I_  *  Til  +  I  has  been  measured  by  Gedeon  et  al.  (1971)  to 

°  *  -14  2 

have  a  cross  section  of  1.05  X  10  cm  . 

Three  body  atomic  recombination  rate  coefficients  typically  lie 
—31  “33  8 

in  the  range  10  to  10  cm  /sec.  We  assume  a  common  value  for 


-32  6 

reactions  6a  and  6b  of  kg  =  1  X  10  ca/sec  and  siailarly  we  assuae 
kjg  =  1  x  IQ"32  ca3/sec. 

2.2.4  Diffusion 

Diffusion  of  Tla  through  a  nuaber  of  buffer  gases  has  been 

measured  by  Aleksandrov  et  al.  (1978)  at  390*C.  The  values  of  D  in  N0 
2  2  o  l 

and  Ar  are  0.47  ca  /sec  and  0.3  ca  / sec,  respectively.  For  our  purposes 

we  take  the  aean  of  these  two  values  and  convert  it  to  the  value 
18  -1  -1 

(DN)  =  4.4  X  10  ca  sec  In  the  absence  of  other  data  we  use  this 

saae  value  for  T1  and  Cl  diffusion, 
o 

In  order  to  evaluate  various  reaction  rates  froa  the  above  cross 
sections  we  aust  assign  a  thermal  velocity  v.  The  assuaption  of  a 
teaperature  of  M700K  is  consistent  with  the  results  of  the  next  section. 
At  this  teaperature  the  species  involved  in  reaction  (2.3)  through  (2.9) 
have  velocities  ranging  froa  2.5  X  10*  ca/sec  to  6.5  x  10*  ca/sec.  For 
siaplicity  we  assign  a  coaaon  value  of  v  =  4  X  10*  ca/sec  to  all  the 
relevant  species. 

In  coaputing  the  diffusion  loss  of  Tla  and  T1q  we  assuae  that 
all  aetastable  are  quenched  at  the  wall  (i.e.  pa  =  1),  while  T1q  is 
taken  to  have  various  values  of  its  sticking  probability  in  order  to 
explore  the  relative  iaportance  of  wall  reactions  in  deteraining  the 
steady  state  density  of  T1q. 

In  the  calculations  presented  in  the  next  section  the 
sensitivity  of  the  results  to  all  the  assuaed  input  paraaeters  is 
examined,  with  the  exception  of  pa  and  v  which  are  fixed  at  the  values 
indicated  above. 

2.8  CALCULATIONS 

The  general  procedure  used  to  solve  for  various  sets  of 
self-consistent  steady  state  conditions  was  derived  in  Section  2.1.  All 
of  the  calculations  reported  here  refer  to  a  right- circular  cylinder 
having  height  hQ  =  20  ca,  and  radius  rQ  =  10  ca.  This  is  believed  to  be 
closed  to  an  optiaua  sise  for  a  single  *aodule,(  possibly  foraing  part 
of  a  larger  systea  in  a  practical  device. 


The  required  concentration  of  Tla  is  specified  by  the  need  to 
have  a  reasonably  high  absorption  probability  for  an  incoaing  535  no 
photon . 

With  a  physical  depth  of  20  ca,  and  an  absorption  cross  section 
-12  2 

of  4  X  10  ca  (C.  H.  Corliss  et  al.  1062),  an  optical  depth  of  unity 

10  -3 

is  achieved  with  [TlaJ  =  1.25  X  10  ca  In  perforaing  various 
exploratory  calculations  we  have  therefore  focussed  attention  on  those 

coabinations  of  input  paraaeters  which  give  values  of  [T1  ]  in  the 

10  3  1 

region  of  1  X  10  /ca  .  Pour  distinct  situations  have  been  considered: 

(1)  Only  T1C1  and  its  products  are  present.  (2)  T1C1  and  its  products 

are  present  with  an  inert  buffer  gas.  (3)  Excess  Clg  is  added  to  the 

T1C1.  (4)  Both  a  buffer  gas  M,  and  excess  Clg  are  added  to  the  aixture. 

In  all  cases  the  teaperature  is  assuaed  to  be  sufficiently  high  for  all 

species  to  be  aaintained  in  the  gas  phase.  This  iaplies  a  lower  liait 

on  the  teaperature,  deterained  either  by  the  prescribed  density  of  T1C1, 

or  by  the  density  of  T1q.  The  saturated  density  of  T1q  varies  froa  1011 

to  10*2  ca  2  for  400*C  <  T  <  450*C.  We  therefore  regard 
12  -3 

[Tl  ]  >  10  ca  as  undesirably  high  and  have  sought  conditions  giving 
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values  closer  to  10  ca 

Exploratory  calculations  deterained  that  a  punping  photon  flux 
17  2 

density  of  1.28  X  10  photons/ca  /sec,  corresponding  to 
2 

100  ailliwatts/ca  of  line  radiation  at  254  na,  was  sufficient  to 

10  -3 

sustain  a  density  of  [Tla]  -  1  X  10  ca  under  appropriate  conditions. 
This  value  of  was  enployed  for  essentially  all  the  calculations 
presented  below.  The  sticking,  or  rather  quenching  probability  for  Tl^ 
striking  the  container  walls  has  been  assuaed  throughout  to  be  unity. 

The  value  of  p^,  the  sticking  probability  of  T1q  was  varied  over  a  wide 
range  to  explore  the  relative  potential  iaportance  of  wall 
reconstitution  of  T1C1.  In  all  cases  the  corresponding  sticking 
probability  of  Cl,  Pj,  *■*  as  a  free  parameter  and  computed  as 

output,  as  explained  in  Section  2.1. 
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2.3.1  No  Buffer  6m  ([M]  ■  0) 

In  this  situation  it  is  appropriate  to  regard  the  density  of 
T1C1  as  the  independent  variable,  and  the  computed  self  consistent 
stoichiometric  densities  are  plotted  vs.  [T1C1]  in  Figure  2.3  for  the 
cases  pj  =  1  and  0.1.  Stoichiometry  is  imposed  on  the  solutions  by 
adjusting  the  input  value  of  Clj  so  as  to  achieve  equal  values  of 
{[T1q]  *  [T1b]}  and  {[Cl]  ♦  SfClg]}-  The  importance  of  reaction  2.3  is 
demonstrated  by  the  two  curves  for  [T1b]  corresponding  to  the  indicated 
values  of  o^,  the  cross  section  for  quenching  of  T1b  by  T1C1.  This 
process  is  dominant  at  high  [T1C1]  densities,  causing  the  density  of  T1 

ID 

to  saturate  at  a  value  where  it  is  quenched  by  T1C1  as  rapidly  as  it  is 

produced  from  it,  i.e.  [T1  ]  -  N  o./y  <?-. 

*  P  .  "  13 

At  sufficiently  low  densities  of  T1C1  (£10  )  the  dominant  loss 

of  T1  is  by  kinetic  free  fall  to  the  container  walls,  where  it  is 
quenched.  Under  these  conditions  the  loss  frequency  depends  only  on  the 
container  sise,  and  on  the  wall  quenching  probability  pa,  set  equal  to 
unity  for  all  present  calculations.  Bence  [T1b]  increases  essential 
linearly  with  [T1C1] ,  due  to  the  linearly  increasing  rate  of  production. 
As  expected,  there  is  no  dependence  of  [T1b]  on  the  T1q  sticking 
probability  p^.  The  [Tl^]  curves  shown  apply  for  both  p^  =  1  and 
Pl  =  0.1. 

Under  the  conditions  of  Figure  2.3,  T1q  is  lost  predominantly  by 
wall  reconstitution,  forming  T1C1.  At  low  values  of  [T1C1] ,  where  "free 
fall"  diffusion  is  occurring,  this  rate  is  strongly  sensitive  to  the 
value  of  p^.  The  resulting  steady  state  density  of  T1q  varies 
approximately  m  (2-p^)/pj,  the  relevant  factor  appearing  in  Eq.  (2.16). 
The  transition  from  "free  fall"  diffusion  to  standard  diffusion  occurs 

9 

in  the  region  where  IX.  ~  A  ,  as  indicated  by  Eq.  (2.15).  For  the 

010  13-3 

present  container  this  occurs  when  N  ~  [T1C1]  -  6  X  10  cm  ,  when 

Pj  =  1,  and  at  an  appropriately  higher  density  when  p^  =  0.1.  This 

transition  causes  the  dependence  of  [T1q]  on  [T1C1]  to  change  from 

linear  to  quadratic.  The  latter  dependence  arises  because  the  diffusion 

loss  rate  of  T1q  is  then  decreasing  inversely  as  N  -  [T1C1] ,  while  its 

production  rate  is  increMing  linearly  with  [T1C1] . 
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Figure  2.3  Dependence  of  the  important  species  densities  on  the  density 

17  2 

of  T1C1,  for  a  fixed  photon  flux  of  1.28X10  photons/cm  /s, 
with  no  added  buffer  gas,  and  with  the  Clj  density  adjusted 
to  give  stoichiometry.  Curves  for  two  values  of  the 
sticking  probability  of  T1q  are  shown,  as  is  the  effect  on 
T1  of  varying  the  cross  section  for  reaction  2.3. 
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The  density  of  Cl  is  not  shorn  explicitly  in  Figure  2.3.  Over 
most  of  the  range  it  is  essentially  the  same  as  [T1q] ,  since 
stoichiometry  requires  [T1q]  +  [Tl^]  =  [Cl]  +  2[Cl2],  and  [Clg]  is 
extremely  small. 

Of  particular  interest  is  the  T1q  density  when  [Tl^]  approaches 

saturation.  This  depends  on  the  T1  sticking  probability.  For  p.  =  1 

11  °-3  1 
the  relevant  T1  density  is  <10  cm  ,  and  with  p.  =0.1  the  value 

12  °  -3  * 

rises  to  ~10  cm  The  latter  may  imply  undesirably  high  temperature, 

and  we  therefore  seek  a  method  other  than  wall  recombination  for 

reducing  the  equilibrium  density  of  T1q.  The  possibility  of 

accelerating  the  three-body  volume  recombination  processes  by  adding  an 

otherwise  inert  buffer  gas  is  explored  in  the  next  subsection. 

2.3.2  Vith  Buffer  Gas,  "M* 

The  effect  of  adding  a  buffer  gas  to  the  mixture  has  been 

investigated  with  a  view  to  enhancing  the  three-body  volume 

recombination  processes  by  which  atomic  thallium  and  chlorine  are 

reconstituted  to  T1C1.  As  an  example  we  show  in  Figure  2.4  the  effects 

on  the  situation  when  a  value  of  =0.1  is  assumed.  The  chain  curves, 

for  no  buffer  gas,  are  reproduced  from  Figure  2.3.  The  addition  of  a 
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buffer  gas  density  of  3  X  10  cm  gives  the  situation  represented  by 
the  full  curves.  As  expected,  the  density  of  Tl^  approaches  saturation 
at  a  much  lower  density  of  T1C1  due  to  the  suppression  of  the  free  fall 
diffusion  loss.  Similarly,  the  densities  of  T1q  and  Cl  increase 
together  in  the  presence  of  the  buffer  gas,  due  to  inhibition  of  their 
transport  to  the  wall  where  their  recombination  is  primarily  occurring. 
There  is  no  evidence  that  the  volume  recombination  processes  have  been 
made  effective  by  this  amount  of  added  buffer  gas,  and  the  overall 
effect  is  that,  in  the  region  where  [Tlo]  is  close  to  its  saturated 
level,  the  density  of  T1q  has  increased  over  the  corresponding  value 
without  the  buffer  gas.  Increasing  the  buffer  gas  density  too  much 
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Figure  2.4  Showing  the  effects  of  adding  a  buffer  gas  and/or  Cl2,  with 

the  T1q  sticking  probability  of  0.1.  The  dot-dash  curves, 

taken  from  Figure  2.3,  correspond  to  no  buffer  gas  and 

stoichiometry.  The  effect  of  adding  01o  to  maintain  its 
12  -3  * 

density  at  1X10  cm  is  shown  by  the  long-dash  curves,  and 

lfi 

the  effect  of  also  adding  3X10  cm  of  buffer  gas  is  shown 

by  the  short  dash  curve.  The  full  curves  correspond  to  a 

16  —3 

buffer  gas  density  of  3X10  cm  with  the  Cl2  density 
adjusted  for  stoichiometry. 


m 


higher  levels  is  a  viable  option  if  pressure  broadening  of  the 

absorption  line  is  needed.  However,  the  use  of  substantially  higher 

buffer  gas  densities  can  cause  noticeable  reduction  of  the  Tl^  density, 

due  to  collisional  quenching.  This  effect  is  seen  in  Figure  2.5,  where 

the  various  densities  are  plotted  as  a  function  of  the  buffer  gas 
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density  for  a  constant  T1C1  density  of  1  X  10  ca  In  this  figure 

the  broken  curves  represent  the  situation  presently  being  considered, 

where  only  T1C1  and  its  stoichiometric  products  are  present  with  the 

buffer  gas.  The  peak  in  the  density  of  T1  occurs  in  the  region  of 
17  -3  0 

[M]  =  6  x  10  ca  ,  where  the  volume  recombination  processes  start  to 
control  the  densities  of  T1q  and  Cl.  At  low  buffer  gas  densities  [T1q] 
is  sensitive  to  its  sticking  probability,  p^,  whose  value  is  indicated 
on  some  of  the  curves  in  Figure  2.5.  Where  not  indicated,  the  curves 
are  insensitive  to  this  parameter.  It  is  clear  from  this  figure  that 
the  only  advantage  of  adding  a  buffer  gas  to  an  otherwise  stoichiometric 
mixture  derived  from  T1C1  is  a  modest  increase  in  the  sustained  density 
of  aetastables,  this  being  achieved  at  the  expense  of  a  significant 
increase  in  the  ground  state  metal  atom  density,  [T1q] . 

2.3.3  With  Bxtra  Chlorine 

As  an  alternative  technique  for  controlling  the  density  of 
atomic  thallium  we  consider  now  the  use  of  excess  chlorine.  To 
investigate  this  possibility  a  series  of  calculations  were  performed  in 
which  the  density  of  Cl^  was  prescribed  at  some  fixed  value,  without 
concern  for  the  overall  stoichiometry.  The  effect  of  adding  a  Cl„ 

12  -3  ^ 

density  of  1  X  10  ca  is  represented  in  Figure  2.4  by  the  long-dashed 

curves,  which  should  be  compared  to  the  dot-dash  curves  representing  the 

situation  with  no  additives.  Note  that  [T1  1  has  decreased  over  the 

o 

whole  range  of  T1C1  densities,  while  [Tl^]  has  increased  at  low 
densities.  As  a  result,  in  the  region  where  [T1  1  now  approaches 

B 

saturation,  there  is  a  substantial  decrease  in  the  T1  density  to  a 

11-3  ° 

value  of  approximately  1  X  10  ca 
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Figure  2.5  Showing  the  dependence  on  buffer  gas  density,  for  a  fixed 

1 3  -3 

T1C1  density  of  1X10  cm  .  The  broken  curves  correspond 
to  stoichiometry  for  Cl9,  while  the  full  curves  are  for  a 
fixed  Cl2  density  of  1X10  cm'  .  Curves  for  two  values  of 
the  T1q  sticking  probability  are  shown  where  there  is  a 
sensitivity  to  this  parameter. 
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The  effect  of  varying  the  density  of  Cl0  over  a  wide  range  is 
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illustrated  in  Figure  2.6  for  a  constant  T1C1  density  of  1  X  10  ca  , 
and  assuming  a  sticking  probability  of  =  0.03  for  T1q  at  the  wall. 

The  broken  curves  represent  the  situation  presently  being  considered, 
where  only  01o  is  added  to  the  T1C1.  We  note  that  densities  of  01o 
greater  than  1  X  10  ca  are  effective  in  reducing  the  density  of  T1q, 

while  the  T1  density  is  essentially  unaffected  until  the  Cl„  density 

B  12  -3  * 

exceeded  1  x  10  cm  .  Thus  it  is  clear  that  the  addition  of  modest 

amounts  of  excess  Clg  provides  control  of  the  steady  state  density  of 

T1  . 
o 

2.3.4  With  Added  Chlorine  and  Buffer  Gas 

A  number  of  calculations  were  performed  to  explore  the  possible 
advantages  of  adding  both  a  buffer  gas  and  excess  chlorine.  In 
Figure  2.4  the  short  dashed  curves  represent  the  situation  where 
[M]  =  3  X  1016  cm  ^  and  [Clg]  =  1  X  10^  ca  ^ .  Compared  to  the  long 
dashed  curves  representing  the  case  with  only  Clg  added,  these  show  a 
minor  increase  in  the  sustained  density  of  T1  accompanied  by  a  similar 

B 

increase  in  [T1  ] .  The  effect  of  varying  the  buffer  gas  density  over  a 

wide  range,  while  maintaining  [Clg]  =  1  X  10  cm  is  shown  by  the  full 

curves  in  Figure  2.5.  The  main  effect  of  the  excess  chlorine  is  that 

its  controls  the  density  of  T1q  over  the  whole  range,  with  only  a  slight 

dependence  of  [T1q]  on  the  sticking  probability  p^  at  the  lower  buffer 

gas  densities.  In  this  region  of  [M]  we  note  that  the  effect  of  the 

excess  chlorine  on  [T1q]  depends  significantly  on  the  magnitude  of  p^. 

When  this  parameter  is  small,  addition  of  Cl0  causes  [T1  ]  to  decrease, 

due  to  the  effect  of  reaction  (2.5).  On  the  other  hand,  if  p^,  is  close 

to  unity,  the  high  atomic  chlorine  density  arising  from  the  added 

chlorine  inhibits  the  transport  of  T1q  to  the  wall,  causing  its 

volumetric  concentration  to  increase. 

The  effects  of  varying  the  chlorine  density  with  a  fixed  buffer 
16  —3 

gas  density  of  3  X  10  ca  is  shown  by  the  full  curves  in  Figure  2.6, 
where  a  small  value  of  p^  =  0.03  has  been  assumed.  In  the  region  of 
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Figure  2.8  Showing  the  effect  of  varying  the  density  of  Cl„  over  a  wide 
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range,  for  a  constant  T1C1  density  of  1X10  ca  ,  and  a 

constant  T1q  sticking  probability  of  0.03.  The  broken 
curves  correspond  to  no  added  buffer  gas,  while  the  full 
curves  are  for  a  fixed  buffer  gas  density  of  3X10  ca  . 
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[Cl0]  -  1  x  10  cm  9  where  the  T1  density  has  been  reduced 
«  o 

substantially  without  significantly  affecting  the  metastable  density,  we 
see  that  additional  buffer  gas  has  very  little  influence  on  either  T1q 

or  T1  . 

a 


2 . S . 6  Sensitivity  Study 

In  the  preceding  subsections  we  have  explored  the  effects  of 
varying  the  physically  controllable  densities,  of  T1C1,  Clg,  and  buffer 
gas  M.  In  addition  the  cross  section  for  reaction  (2.3),  and  the  T1q 
sticking  probability  Pj  have  been  varied.  It  appears  from  these  results 
that  there  exists  a  rather  extensive  parameter  space  within  which  it  is 
possible  to  sustain  the  desired  density  of  metastables  with  modest 
optical  pumping  power  levels.  It  is  important  to  establish  that  this 
conclusion  does  not  rest  on  some  fortuitous  circumstance  dependent  on 
the  assumptions  made  regarding  the  many  other  parameters  involved  in  the 
model.  To  this  end,  additional  calculations  have  been  performed  where, 
in  turn,  each  input  parameter  was  first  increased  a  factor  of  ten,  and 
then  decreased  a  factor  of  ten  from  its  baseline  value.  The  chosen 
baseline  case  parameters  follow  from  the  considerations  discussed  in 
Section  3,  or  from  the  results  presented  above. 

The  results  of  this  sensitivity  study  are  presented  in 

Table  2.2.  The  parameters  involved  are  listed  in  the  first  column, 

followed  in  the  second  column  by  the  values  assigned.  The  final  three 

columns  give  the  atomic  densities  of  interest.  The  first  row 

corresponds  to  our  baseline  case  around  which  the  variations  are  being 

explored.  For  this  case  each  parameter  listed  has  a  value  equal  to  the 

geometric  mean  of  the  two  values  listed  in  the  table.  For  example,  for 

this  case,  ff.  =  3.5  x  10  cm, as  indicated  in  the  first  row,  while 
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ff2  =  1.5  x  10  cm  etc. 

Examination  of  this  table  reveals  certain  anticipated 
sensitivities.  The  sustained  density  of  metastables  is  more  or  less 
proportional  to  the  parameters  and  N^,  which  control  the  rate  of 
reaction  (2.1).  Reaction  (2.2)  plays  a  significant  role  in  determining 


the  density  of  T1q,  but  has  essentially  no  effect  on  the  aetastable 
density.  The  converse  is  true  of  reaction  (2.3),  whose  effect  was 
explored  somewhat  in  Figure  2.3.  Reaction  (2.4)  is  quenching  of 


metastables  by  Clg,  and  the  effects  of  varying  a ^  are  consistent  with 
the  effect  of  varying  [Clg]  illustrated  in  Figure  2.6.  Reaction  (2.5) 
plays  a  dominant  role  in  controlling  the  density  of  T1q,  whose  density 
is  seen  to  be  essentially  inversely  proportional  to  a g.  The  only  other 


significant  sensitivity  is  to  the  pumping  photon  flux  density,  N^. 
Changes  in  the  magnitude  of  this  parameter  cause  changes  in  the 


densities  of  T1q  and  Tl&  which  are  similar  to  the  effects  of  the  changes 


made  in  a ^  and  respectively.  Since  is  subject  to  external  control 
it  can  in  principle  be  adjusted  to  compensate  for  any  error  in  the 


assumed  value  of  a j. 


Table  2.2 


[Cl] 

[T1J 

-3 

-3 

-3 

ca 

ca 

ca 

Value 

1B+11 

1B+15 

1B+09 

3. SB-17 

1.58 

0.46 

9.25 

3. SB-18 

0.58 

0.28 

0.04 

3. SB-16 

11.53 

1.22 

86.82 

1. SB-18 

1.15 

0.30 

0.30 

1. SB-16 

5.84 

0.87 

8.04 

l.OE-15 

1.58 

0.46 

34.54 

1. OB-13 

1.58 

0.46 

1.11 

1. SB-15 

1.58 

0.46 

10.39 

1. SB-13 

1.58 

0.46 

4.41 

1. OB-15 

12.50 

0.41 

9.25 

1. OB-13 

0.16 

0.46 

0.25 

T. OB-33 

1.58 

1.42 

8.58 

1. OB-31 

1.58 

0.14 

0.40 

5. OB-17 

1.58 

0.46 

0.26 

5. OB-15 

1.58 

0.46 

0.20 

1 .08-20 

1.58 

0.46 

6.46 

1. OB-18 

1.58 

0.46 

7.50 

1. OB-18 

1.58 

0.46 

0.57 

1. OB-16 

1.58 

0.46 

6.04 

1. OB-33 

1.58 

0.46 

6.25 

1.0E-31 

1.57 

0.45 

0.25 

2.2E+18 

1.60 

0.46 

0.26 

8.8E+18 

1.54 

0.45 

9.24 

1.0E+00 

1.58 

0.46 

6.25 

1.0E-02 

1.59 

0.46 

0.25 

1.3E+16 

0.16 

0.14 

0.05 

1.3E+18 

15.76 

1.42 

85.36 

3.  THE  Tl-Cs  SYSTEM 

3  COLLISIONAL  TRANSFER  BXCITATION 

The  collisional  transfer  excitation  of  the  setastable  state  of 
thallium  atoms  involves  introducing  into  thalliua  vapor  another  species 
with  resonance  states  in  the  vicinity  of  the  setastable  state  level  of 
Tl,  to  which  the  collisional  transfer  from  those  resonance  states  can 
easily  take  place.  Figure  3.1  shows  a  partial  energy  level  diagram 
which  includes  the  thallium  states  involved  as  well  as  some  appropriate 
cesium  states.  The  radiation-trapped  Cs(6^P)  states  have  long  lifetimes 
and  then  can  transfer  energy  to  ground  state  Tl  atoms  to  excite  them  to 
the  metastable  state  Tl (^*3/3)  collisions.  The  transfer 

excitation  rate  is  a  strong  function  of  the  energy  defect  between  the 
two  states  involved  (Csajkowski  et  al. ,  1073),  in  this  case  about 
0.4  eV.  The  excitation  transfer  cross  sections  are  estimated  to  be 
2  x  10  ^  cm^  for  the  Tl-Cs  system. 

Table  3.1  lists  the  various  reactions  and  their  cross  sections 
relevant  to  the  performance  of  an  ALF  employing  Cs-Tl  mixtures.  Those 
cross  sections  not  available  in  the  literature  have  been  estimated.  The 
sustained  density  of  Tl^  depends  primarily  on  the  balance  between 
production  by  Reaction  (3.2)  and  loss  by  reactions  (3.5)  and  (3.6).  The 
corresponding  density  of  Cs*(8P)  that  needs  to  be  sustained  is  thus 
given  by 

.  Cs  oA 

[Cs*(8P)]  =  Tl  -I— 2,— *  ♦  -2  (3.0) 

T1„  2 

OJ  it 

The  pumping  power  required  to  sustain  this  density  is  determined 
primarily  by  Reaction  8  and  the  radiative  loss  of  Cs*(6P)  modified  by 


Table  3.1 


List  of  Reactions  and  Cross  Sections 


Reactions 


Cross  Sections  (c«  ) 


(1)  Cs(6S)  +  hv(890)  ♦1  Cs*(6P) 


0^  “  1.7  x  10' 


2  _,b/2 


(2)  Cs  (6P)  +  Tl(  p1/2)  *  T1  (  P 3/2)  +  Cs(6S)  +  AE  Oj  ■  2  x  10' 


(3)  T1"(2P1/,)  +  h  v(535)  >3  T1*(7S) 


-  4.9  x  10' 


(4)  T1*(7S)  «■*  Tl(2Pl/2)  +  h v(378) 


3  x  10“12 


5 


(5)  Tlm(^P3/2)  +  Cs(6S)  ♦  Tl(£P1/2)  +  Cs(6S)  +  AE  a,  -  5  x  10' 


(6)  Tl“(zP1/9)  +  Tl(‘p.  ,,)  2T1(4P, /-)  +  AE 


^  1T1 t 2. 


o6  ■  5  x  10' 


(7)  Cs  (6P)  +  Tl(2P1/2)  *  Cs(6S)  +  Tl(2 7P1/2)  +  AE  o7  -  5  x  10" 


(8)  Cs*(6P)  +  Cs(6S)  2Cs(6S)  +  AE 


(Jg  ■  2.  1  X  10 


(9)  Cs*(6P)  +  Cs*(6P)  *9  Cs**  (6D)  +  Cs(6S)  +  AE 


1.5  x  lO"13  -  2  x  10' 


(10)  Cs*(6P)  +  h y( 37 7 )  -10  Cs+  +  e 


(11)  Cs*(6P)  +  T1*(7S)  *nCs+  +  Tl(2P1/2)  +  e 


<jj  j  •  5  x  10 


(12)  Cs*(6P)  +  Cs*(6P)  ■►I2Cs2+  +  e 


1.5  x  10' 


1 1 


K'y 


b 


sr 


V> 


& 
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resonance  radiation  trapping.  Together,  these  processes  impose  an 
effective  lifetime  of  on  the  Cs*(6P)  atoas.  This  effect  has  been 

studied  in  pure  Cs  by  recording  the  growth  and  decay  of  the  fluorescence 
signals  from  both  resonance  states  following  the  pumping  of  one  or  other 
of  the  states  by  a  short  pulse  fron  a  tunable  dye  laser. 

S.2  TIB  BXPBKIMBNT 

A  long  cylindrical  shape  was  chosen  for  the  experinental  cell  to 

approximate  the  atonic  line  filter  geometry.  For  the  apparent  lifetime 

study  of  the  CsCSPjyg)  u<^  Cs^Pg/g)  states  we  have  used  for  the 

excitation  a  laser  pulse  which  is  short  compared  to  the  natural 

lifetimes.  This  simplifies  the  analysis  of  the  measured  time  dependent 

fluorescence  of  the  parent  state  and  the  excitation  transferred  state. 

The  cesium  vapor  density  covered  in  the  present  experiment  ranges  from 
11-3 

10  cm  ,  where  the  Doppler  effect  dominates  the  line  shape,  to 
14-3 

10  cm  ,  where  collisional  effects  dominate  the  line  shape.  Ve 
describe  first  the  physical  characteristics  and  the  preparation 
procedure  of  the  cesium  cell,  and  we  then  describe  the  experimental 
setup  and  the  method  by  which  the  data  are  collected  and  prepared  for 
analysis . 

3.2.1  The  Cesium  Cell 

The  cesium  cell  is  a  cylinder  made  of  304  stainless  steel,  23  cm 
long  and  inside  diameter  11.5  cm,  simulating  one  module  of  the  atomic 
line  filter.  It  is  shown  schematically  in  Figures  3.2  and  3.3.  It  has 
three  optical  ports,  each  made  of  M3  mm  thick  sapphire  brased  to  a  kovar 
sleeve  and  mounted  on  a  standard  size  Varian  vacuum  flange.  The  two  end 
ports  are  ~25  mm  in  diameter,  allowing  the  excitation  laser  pulse  to 
pass  through  the  cesium  cell.  The  third  port,  M50  mm  in  diameter,  is 
mounted  approximately  mid-way  on  the  side  wall  and  is  oriented  90*  to 
the  first  two.  An  ~2  ft*  thick  layer  of  nickel  is  coated 
electrolytically  on  the  sapphire  -aetal  joints  to  prevent  their 
deterioration  from  cesium  attack.  The  copper  gaskets  for  sealing  the 
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sapphire  window  flange  to  the  cell  body  are  also  plated  with  ~200  A  of 
nickel. 

The  arrangement  used  to  vacuum  process  the  cell  is  shown  in 

Figure  3.2.  A  13  mm  diameter  U-tube  to  be  used  as  the  cesium  reservoir 

is  welded  to  the  center  bottom  of  the  cesium  cell.  The  other  end  of  the 

U-tube  has  a  15  cm  long  straight  section  of  thin  wall  (~0.25  am) 

stainless  steel  tubing  connecting  the  cesium  cell  to  a  conventional  high 

vacuum  system.  While  pumping,  the  cell  is  baked  at  ~400*C  for  12  hours 

-8 

or  more  until  a  vacuum  of  M10  torr  is  attained.  Prepurified  cesium  in 
a  glass  capsule  is  then  opened  and  about  2  grams  of  it  is  distilled  into 
the  cell.  The  thin  wall  stainless  steel  tubing  is  pinched  to  form  a 
sealed  bridge,  and  a  controlled  amount  of  current  is  passed  through  the 
bridge  to  fora  a  resistive  weld,  thereby  capping  off  and  separating  the 
cesium  reservoir  from  the  vacuum  system.  The  cesium  cell  is  then 
transferred  to  the  experimental  station  described  in  the  next 
subsection. 

3.2.2  Experimental  Apparatus 

The  apparatus  used  to  determine  the  effective  lifetimes  of  the 
CsCBPjyg)  and  s^a^es  is  shown  schematically  in  Figure  3.3. 

There  are  three  major  components:  (1)  the  cesium  cell  which  provides 
the  proper  container  shape  and  the  necessary  cesium  vapor  density, 

(2)  the  tunable  dye  laser  that  provides  pulse  excitation  of  the  cesium 
vapor,  and  (3)  the  components  for  detection  and  display  of  the  time 
dependence  of  the  direct  and  the  excitation  exchange  fluorescence 
signals. 

After  processing,  the  cesium  cell  is  housed  in  a  specially 
constructed  oven  made  of  8.4  cm  thick  fire  bricks,  the  inside  lining  of 
which  is  made  of  1.5  mm  aluminum  sheets.  Three  sets  of  heaters  are 
used,  one  set  on  the  ceiling  panel,  one  set  on  the  bottom  panel  and  one 
set  on  the  side  walls.  Ten  thermocouples  strategically  attached  to  the 
cesium  cell  body  monitor  its  temperature  distribution.  Heater  powers 
are  adjusted  such  that  the  cell  body  temperature  is  always  higher  than 
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Figure  3. 


L  =  Lens 
SF  =  Step  Filter 
F  =  Narrow  Band  Filter 
ND  =  Neutral  Density  Filter 
PMT  =  Photomultiplier 


3  The  optical  arrangement  used  to  pump  and  view  the  Cs 
fluorescence. 


that  of  the  cesium  reservoir,  the  difference  being  a  few  degrees  at  low 
cesium  densities,  to  as  high  as  30*C  at  high  densities.  Once  the  cell 
body  reaches  a  predetermined  temperature,  a  temperature  controller 
regulates  the  power  to  the  heaters.  Temperature  stability  of  the  cesium 
cell  can  normally  be  maintained  to  within  1*C. 

There  are  three  double-walled  quarts  windows  5  cm  in  diameter  in 
the  oven’s  side  walls,  aligned  with  the  cesium  cell’s  sapphire  windows. 
Two  of  these  windows  are  for  the  entrance  and  exit  of  the  excitation 
laser  pulse,  and  the  third  one  90*  to  the  other  two  is  for  the 
observation  of  the  fluorescence. 

The  cesium  vapor  is  optically  pumped  by  a  pulsed  and  spectrally 
tunable  dye  laser  (Quanta  Ray  Model  DCR/PDL)  pumped  by  a  Nd:YAG  laser 
which  is  flash-lamp-pumped,  Q-switched  and  frequency  doubled  to  532  nm 
for  the  dye  laser  pumping.  With  Exciton  LDS  867  dye  and  at  10  Hz  the 
dye  laser  output  is  ~200  nJ  per  pulse  at  wavelengths  852  nm  or  894  nm. 
The  laser  pulse  duration  is  ~9  ns  and  the  spectral  width  is  ~0.4  cm*. 
The  laser  beam  axis  is  aligned  collinearly  with  that  of  the  cesium  cell. 
The  unstable  oscillator  resonator  design  of  the  Nd:YAG  laser  results  in 
a  donut  shaped  pump  beam.  This  in  turn  produces  a  donut  shaped  dye 
laser  beam  measuring  ~8  mm  outside  diameter  and  ~3  mm  inside  diameter. 
Because  the  dye  fluoresces  strongly  at  the  pump  frequency  (532  nm) ,  a 
Schott  RG  830  red  filter  is  placed  at  the  entrance  port  of  the  cesium 
cell  to  preclude  any  possible  interference  with  the  time  dependent 
resonance  photon  transport  processes.  The  dye  laser  output  is  usually 
measured  with  a  Scientech  362  power  meter  at  a  point  between  the  RG  830 
red  filter  and  the  entrance  port  to  the  cesium  cell,  prior  to  and 
shortly  after  each  experimental  run.  When  tuned  and  stabilized  properly 
we  find  no  significant  change  of  laser  output  at  the  set  wavelength  over 
a  period  of  an  hour.  The  laser  excitation  pulse  energy  was  reduced  by 
predetermined  amounts  with  a  calibrated  neutral  density  filter  to  check 
for  any  dependence  of  the  time  dependent  resonance  fluorescence  on  power 
density.  For  the  laser  excitation  pulse  energy  used  in  the  present 
experiment  there  is  no  noticeable  effect  on  the  measured  effective 
lifetime  of  either  the  Cs^Pjyg)  or  the  CstbPg/g)  state. 


3.3  852  na  AND  894  na  FLU01ESCBNCB  FROM  PULSB  EICITED  Cs(0P1/9)  AND 

Cs (®Pgy2>  STATES 

3.31  Bxperiaental  Conditions 

After  reaching  a  prescribed  cesitia  vapor  density  in  the  cell, 
the  dye  laser  system  was  turned  on  and  tuned  to  either  852  nm  or  894  nm 
wavelength  for  resonance  excitation  of  the  CstBPg^g)  or  the  CstBPj^g) 
state.  To  avoid  initial  nonuniform  excitation  along  the  axial  direction 
of  the  long  cesium  cell  we  set  the  dye  laser  pulse  energy  at  a  level 
exceeding  the  burn-through  or  column  saturation  value.  The  latter  is 
defined  as  that  power  density  which  provides  approximately  one  resonance 
photon  for  each  atom  in  the  laser  beam  column  per  radiative  lifetime 
[Huennekens  and  Gallagher  1983].  That  is, 


col . sat. 


“•  h"l2l  [“Wr 


where  N  is  the  cesium  atom  density,  L  is  the  column  length,  hi/^  is  the 
resonance  photon  energy,  is  the  natural  lifetime  of  the  upper  state 

[CsCBP^g)  or  Cs^Pjyg)]  >  and  ^laser  ^atom  are  ^®  l48®1*  line 

width  and  atomic  absorption  line  width,  respectively.  The  ratio 

tv  .  /tv,  represents  the  fraction  of  laser  photons  that  can  be 
atom'  laser  r  r 

absorbed  by  the  atoms  and  is  set  equal  to  1  if  Aj/  .  i  tv,  .  If  the 
'  atom  laser 

laser  pulse  duration  is  shorter  than  the  atomic  natural  lifetime, 

should  be  replaced  by  rT  in  Equation  (3.13).  For  example  the  absorption 
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equivalent  width  for  852  nm  at  N  =  6  X  10  cm  is  approximately 
8.2  GHz  [see  next  section]  while  the  laser  linewidth  is  ~12  GHz.  The 
laser  pulse  duration  of  9  nsec  is  much  smaller  than  the  natural  lifetime 
of  30.5  nsec  for  the  Cs^Pg/g)  state.  Therefore 


I  .  .  =  2.6  X  103  Watt/cm2 

column  sat 


This  is  about  a  factor  a  fifty  smaller  than  the  laser  pulse  intensity  of 
1.1  x  10  fatt/ca  at  200  mJ  per  pulse  when  unattenuated.  Therefore  for 
all  practical  purposes  we  have  a  column  saturation  condition  (or 
burn-through  condition) ,  the  cesium  atoms  within  the  cell  are  more  or 
less  uniformly  excited  in  the  axial  direction  along  the  laser  beam  path. 

To  facilitate  the  data  reduction  it  is  desirable  that  the 
spatial  distribution  of  the  resonance  excitation  relax  quickly  to  the 
fundamental  mode  [Holstein  1047,  1951].  This  mode  is  strongly  excited 
initially  by  the  on-axis  pumping,  and  higher  order  modes  will  in  any 
event  decay  much  more  rapidly  [van  Trigt  1960] .  The  assumption  that 
only  the  fundamental  mode  contributes  to  the  measured  decay  rates  is 
therefore  reasonable,  and  is  consistent  with  the  observed  double 
exponential  character  of  the  decays. 

8.3.2  Data  Acquisition 

After  column  saturation  pulse  excitation  of  the  cesium  vapor  to 
the  (6^x/2^^3/2^  8tmte  with  894/852  nm  laser  light  pulse,  the  time 
dependence  of  the  direct  894/852  nm  fluorescence  and  of  the  excitation 
transferred  852/894  nm  fluorescence  are  measured.  A  75  mm  diameter, 

130  mm  focal  length  lens  is  placed  close  to  the  double-walled  quartz 
oven  window  to  collect  the  fluorescent  light  from  the  mid-section  of  the 
cell  and  focus  it  on  the  SI  photocathode  of  the  K2276  photomultiplier 
(PUT) .  Bach  line  is  isolated  by  a  narrowband  (M10  nm  FWHM)  interference 
filter  having  a  peak  transmittance  of  ”*45%.  Linearity  of  the  PUT 
response  is  checked  with  a  set  of  calibrated  neutral  density  filters, 
and  is  found  to  be  adequate.  The  electron  transit  dispersion  time  of 
the  PUT  is  reported  to  be  a  few  nanoseconds  which  is  adequate  for  the 
present  purposes.  Output  of  the  PUT  is  fed  to  a  Tektronix  7834  storage 
oscilloscope  equipped  with  a  7A26  dual  trace  amplifier  that  is 
terminated  with  a  100  Ohm  resistor.  The  bandwidth  of  the  7A26  amplifier 
is  dc  to  200  UHz,  when  terminated  with  a  50  Ohm  load.  The  recorded 
traces  are  photographed.  Subsequently  the  photographs  were  enlarged  to 
a  size  convenient  for  digitizing  using  a  digitally  interfaced  plotter 
(Tektronix  Uodel  No.  4662). 
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Figures  3.4,  3.5  and  3.6  are  typical  oscilloscope  traces  of  the 
direct  and  excitation  transfer  fluorescence  at  three  different  cesium 


densities:  2.55  X  10**  cm  3,  7.3  X  10*2  cm  3,  and  2.3  X  10**  cm  3, 


respectively.  The  dye  laser  pulse  width  is  ~8  ns  and  the  laser  pulse 
repetition  rate  is  10  Hs.  Subsequent  analysis  shows  that  these  three 
examples  represent  different  regimes.  At  the  lowest  cesium  density  the 
decay  is  dominated  by  radiation  trapping  of  the  Doppler  broadened  lines, 
at  the  intermediate  density  the  line  is  collisionally  broadened,  while 
at  the  highest  density  collisional  self  quenching  dominates  the  observed 
decays . 


The  energy  difference  between  the  6P3^2  U(^  the  6P^2  state  is 


only  0.069  eV.  Therefore  the  following  exchange  collisions  are 
effective  in  the  temperature  range  of  our  experiment,  318*K  to  457*K, 


Cs(8P3^2)  +  Cs(6Sjy2)  ♦  CsCeP^)  +  Cs^S^j)  ♦  AE 


(3.14) 


where  AE  is  the  energy  defect,  and  can  be  taken  up  by  the  thermal  energy 
of  the  participants.  Therefore  experimentally  the  intensity  of  the 
excitation  exchange  fluorescence  should  be  a  linear  function  of  the 
direct  fluorescent.  This  is  shown  in  Figure  3.7  where  the  dye  laser 


wavelength  is  tuned  to  the  6Pjy2  state  (894  nm) ,  and  the  direct  (894  nm) 


and  the  excitation  exchange  (852  nm)  fluorescence  are  both  measured  as  a 


function  laser  power  keeping  the  cesium  vapor  density  at 


7.3  X  1012  cm-3. 


Figures  3.8  and  3.9  show  examples  of  the  late  time  behavior 


12  -3 

observed  at  a  density  of  7.3  X  10  cm  following  excitation  with  a 


laser  pulse  power  of  ~2.2  X  10*  Watts,  or  200  pJ  per  pulse  at  852  nm  and 
884  nm  respectively.  In  these  cases  simple  linear  fits  to  the  data  are 
indicated.  As  predicted  by  the  analysis  of  the  next  section,  all  have 
essentially  the  same  time  constant,  in  this  case  27  pa.  We  have  checked 
for  any  effects  of  the  laser  pulse  power  on  the  fluorescent  time 
constants.  Figure  3.10  is  a  plot  of  the  late  time  fluorescent  decays  of 
894  nm  light  at  200  mJ/pulse  and  63  mJ/pulse  excitation  levels.  Here 
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Figure  3. 


4  Oscilloscope  traces  of  the  direct  and  excitation  transfer 

fluorescence  signals  observed  at  a  cell  temperature  of  332  K 

and  a  reservoir  temperature  of  317  K,  giving  a  Cs  density  of 
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2.55X10  cm  .  The  pump/fluorescence  wavelengths  are 

(a)  894  nm/894  nm,  (b)  894  nm/852  nm,  (c)  852  nm/852  nm,  and 

(d)  852  nm/894  nm. 
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Figure  3.6  Oscilloscope  traces  of  the  fluorescence  signals  observed  at 
a  cell  temperature  of  457  K  and  a  reservoir  temperature  of 
427  K,  giving  a  Cs  density  of  2.3X10^  cm  ^ .  The 
pump/f luorescence  wavelengths  are  (a)  894  nm/894  nm, 

(b)  894  nm/852  nm,  (c)  852  nm/852  nm,  and  (d)  852  nm/894  nm 
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Figure  3. 


Direct  Fluorescent  Intensity  ( 894  nm) ,  Relative 


7  Correlation  of  direct  (894  nm)  and  excitation  transferred 
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(852  nm)  fluorescence  at  cesium  density  of  7.3X10  cnt  \ 
variable  dye  laser  pulse  energies. 
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Figure  3.8  Seai-log  plot  of  the  direct  (852  na)  and  the  excitation 

transferred  (894  na)  fluorescence  intensity  versus  tiae  at 

late  tiaes  after  pulse  excitation  of  the  cesiua  vapor.  The 
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cesiua  vapor  density  is  7.3X10  ca 
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and  in  general  we  see  no  significant  dependence  of  the  decay  time 
constant  on  laser  power. 

The  next  section  is  devoted  to  the  detailed  analysis  of  the  time 
dependent  fluorescence  from  both  the  ®^3/2  s^a^es>  taking  into 

account  the  effects  of  hyperfine  splitting,  line  shape,  excitation 
exchange,  and  self -quenching. 

3.4  DATA  ANALYSIS 

In  this  section  we  describe  the  steps  involved  in  extracting  the 
important  basic  parameters  from  the  time  dependent  fluorescence 
measurements  described  in  the  previous  section.  We  first  present  the 
general  analytical  solution  for  the  three  level  system  involved.  The 
approach  used  to  optimise  the  fitting  of  the  double  exponential 
functions  to  the  data  is  then  described.  The  dependence  of  the 
resulting  two  time  constants  on  the  cesium  density  is  then  fitted  to  a 
model  which  takes  account  of  the  collisional  quenching  of  the  resonance 
states  and  of  trapping  of  the  resonance  radiation.  The  latter  involves 
the  transition  from  a  predominantly  Doppler  broadened  line  to  a 
predominantly  collision  broadened  line,  the  treatment  of  which  is 
further  complicated  by  the  presence  of  hyperfine  components  whose 
spacing  is  comparable  to  the  line  widths,  depending  on  the  Cesium 
density.  The  method  used  to  represent  this  situation  is  described  in  a 
separate  subsection.  In  the  final  subsection  theory  and  experiment  are 
reconciled  by  choosing  appropriate  values  for  the  collisional  energy 
exchange  and  quenching  cross  sections,  which  constitute  the  "unknowns' 
in  the  present  study. 

3.4.1  Time  Dependent  Collisional  Mixing  Model 

The  situation  to  be  analyzed  is  represented  in  Figure  3.11.  The 
analysis  has  been  presented  by  a  number  of  previous  authors  [Rosser 
et  al.  (1069);  Stephenson  and  Bradley  Moore  (1972);  Yardley  (1980)]. 
Therefore  we  give  only  an  outline  of  the  method.  The  notation  employed 
here  and  elsewhere  in  this  section  is  summarized  in  Table  3.2.  The 
continuity  equations  for  the  populations  of  levels  1  and  2  are 
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dNj/dt  =  -N^Kj  +  K12)  +  N2K21 
dN2/dt  =  -N2(K2  *  K^)  ♦  NlKl2 


(3.15) 


These  equations  are  solved  straightforwardly  by  further 
differentiation  and  elimination  of  one  or  other  of  the  variables.  The 
resulting  second  order  differential  equation  is  then  amenable  to 
standard  methods  of  solution.  In  the  present  context  the  solutions  of 
interest  satisfy  the  boundary  condition  that,  at  t=0,  only  one  upper 
level  has  been  populated.  Designating  this  to  be  level  1,  the 
subsequent  behavior  is  given  by 


where 


NjM/N^O)  =  Af  exp(-t/rf)  +  As  exp(-t/rg) 
N2(t)/N1(°)  =  A  [exp(-t/Tg)  -  exp(-t/T{)] 


*  -  >W<V1  -  T»'1) 


(3.16) 


(3.17) 


(3.18) 


*£=(*!*  K12  -  V1)/^,-1  -  r/1) 


(3.19) 


*=  -  (V1  -  Ki  -  -  V‘> 


(3.20) 


r,’1  =  1  {(Kj  .  K12  .  .  K2)  ♦  ((Kj  .  K12  -  K21  -  K2); 


*  «12*21 


,1/2} 


(3.21) 


V1  -  \  {<*1  *  *12  *  K21  *  V  *  *12  -  *21  -  V2  *  4K12K2ll1/2) 


(3.22) 
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Table  3.2  Listing  the  symbols  used  throughout  Section  3  of  this 
report,  together  with  their  definition,  and  the 
equation  where  each  one  is  introduced  or  defined. 

Symbol 

Description 

Eqn. 

a 

Voigt  parameter,  proportional  to  ratio  of  the  collision 
broadened  width  to  the  Doppler  width. 

3.30 

► 

Cfl 

Amplitude  of  the  fast,  slow  components  of  the  double 
exponential  decay. 

3.16 

f 

Actual  number  of  hyperfine  components  of  the  line 

f 

e 

Effective  number  of  hyperfine  components. 

3.27 

F 

Total  angular  momentum  quantum  number  of  a  hf  level 
including  nuclear  spin. 

3.46 

g 

Escape  factor. 

3.24 

*c 

Escape  factor  in  the  presence  of  collision  broadening 
alone . 

3.25 

«d 

Escape  factor  in  the  presence  of  Doppler  broadening 
alone . 

3.25 

gcd 

Escape  factor  associated  with  Doppler  absorption, 
collision  emission. 

3.25 

gu,l 

Degeneracy  of  the  upper,  lower  state. 

3.31 

I 

Nuclear  angular  momentum  quantum  number. 

3.46 

J 

Total  angular  momentum  quantum  number  of  orbital 
electrons  including  electron  spin. 

3.46 

k 

Boltzmann  constant. 

3.30 

k 

o 

Absorption  coefficient  at  line  center,  calculated 
assuming  no  hfs. 

3.31 

k. 

l 

Absorption  coefficient  at  line  center  for  the  i’th 
hyperfine  component. 

3.33 

Kl,2 

Net  loss  coefficient  of  level  1,2  to  the  ground  state. 

3.15 

K12 

Loss  coefficient  of  level  1  to  level  2. 

3.15 

K21 

Loss  coefficient  of  level  2  to  level  1. 

3.15 

40 


Mass  of  the  Cs  atom. 


Total  number  density  of  Cs  atoms. 

Number  density  of  Cs  atoms  in  excited  state  1,2. 

Radius  of  the  cylindrical  container. 

Time  at  which  ^(t)  maximizes. 

Reservoir  temperature,  K. 

Cell  temperature,  K. 

Transmission  coefficient  for  optical  depth  p. 

Most  probable  atomic  speed,  (2kTc/M) 

Average  center-of-mass  thermal  velocity. 

Fractional  relative  intensity  of  the  i’th  hyperfine 

Separation  of  adjacent  hyperfine  components. 

Wavelength  at  line  center. 

Optical  frequency  at  line  center. 

Quenching  cross  section  of  level  1,2  to  the  ground 
state. 

Cross  section  for  excitation  transfer  from  level  1  to  2. 

Effective  lifetime  due  to  collision  broadening. 

Time  constant  of  the  fast,  slow  components  of  the  double 
exponential  decay. 

Natural  radiative  lifetime  of  the  excited  level. 

Displacement  from  line  center  where  k.R  =  1. 

(Absorption  half-width) . 

Absorption  half-width  of  a  collision  broadened  line. 
Absorption  half-width  of  a  Doppler  broadened  line. 


SO 


(3.23) 


We  note  from  Eq.  (3.17)  that  Nj(t)  peaks  at  time 


-1 


-U 


t_  =  ln(rB/rf)  /  (rf  -  rfl  *) 


In  the  above  the  subscripts  *f*  and  *s(  indicate  the  fast  and 
slow  components  of  the  double  exponential  behavior. 

The  solutions  corresponding  to  pumping  level  2  are  obtained  from 
the  above  by  simply  interchanging  1  and  2  in  the  various  subscripts. 

This  has  no  effect  on  and  Tg,  and  hence  tp  is  also  unchanged.  The 
ratio  Af/Ag,  however,  becomes  equal  to  the  reciprocal  of  its  former 
value.  We  conclude  that,  irrespective  of  the  level  originally  pumped, 
the  time  dependent  fluorescence  of  both  levels  1  and  2  can  yield  three 
primary  quantities:  r^,  Tg  and  the  ratio  Aj/Ag .  In  practice  it  is 
possible  to  extract  the  ratio  A^/A g  with  useful  accuracy  only  when  its 
value  is  reasonably  close  to  unity,  and  and  Tg  are  substantially 
different. 


3.4.2  Curve  Fitting  to  N^(t)  and  Nj(t) 

The  procedure  adopted  to  extract  the  values  of  and  Tg  from 
the  digitized  measured  fluorescence  measurements  is  simpler  to  apply  in 
the  case  of  the  Ng(t)  data,  since  in  this  case  the  data  and  the  fitted 
curve  can  both  be  normalized  to  the  peak  value,  and  the  problem  reduces 
to  one  of  simply  fitting  the  shape.  This  was  done  with  an  interactive 
screen  plotting  routine  using  commercially  available  software  on  a 
personal  computer.  Having  digitised  the  data,  it  consisted  of  a 
sequence  of  points.  For  each  such  point  a  corresponding  fitted  point  is 
calculated  using  the  current  best  estimate- of  and  T“,  With  the 

resulting  curves  being  displayed  automatically  in  both  linear  and 
logarithmic  form.  The  program  also  tracks  the  standard  deviation 
between  the  data  and  the  fit,  calculated  in  terms  of  their  logarithms. 

In  seeking  the  best  fit  the  values  of  r ^  and  Tg  were  allowed  to  assume 
any  values  satisfying  the  constraint  <  Tg.  The  convergence  to  the 


fix 


best  values  was  achieved  with  the  help  of  a  routine  which  allowed  one  or 
other,  or  both  of  these  to  be  varied  in  a  prescribed  sequence,  and 
halted  at  the  local  minimum  in  tbe  standard  deviation  of  the  data  from 
the  fitted  curve.  In  the  final  optimization  small  adjustments  were 
allowed  in  the  zero  of  the  time  scale  attached  to  the  data  by  the 
digitization  procedure. 

The  same  general  procedure  was  applied  to  the  fitting  of  N^(t) 
data  except  that  in  this  case  it  was  also  necessary  to  fit  the  amplitude 
normalization. 

Examples  of  both  types  of  fitted  data  are  shown  in  Figure  3.12. 
The  values  of  and  Tg  used  to  compute  each  fitted  curve  are  indicated 
in  the  figure.  Also  indicated  is  the  quality  of  the  fit  as  measured  by 
the  standard  deviation,  quoted  here  as  a  percentage.  In  the  examples  of 
Nj(t)  type  data  the  ratio  A^/Ag  given  by  the  fit  is  also  indicated  on 
the  figure.  Note  that  the  four  curves  shown  in  Figure  3.12  constitute  a 
set,  in  that  all  four  correspond  to  the  same  set  of  experimental 
conditions  of  reservoir  and  cell  temperatures.  Thus  in  principle  they 
should  all  give  the  same  values  for  and  T  ,  and  the  product  of  the 
two  values  of  the  ratio  Aj/Ag  should  be  unity.  For  purposes  of  further 
analysis  we  simply  calculate  the  mean  values  of  and  Tg  and  regard 
these  as  representative  of  the  measurements  at  these  particular 
conditions. 

The  various  sets  of  data  corresponding  to  other  temperatures 

were  treated  in  the  same  way,  giving  in  all  18  pairs  of  values  of  r,  and 

11  * 

r  corresponding  to  Cesium  densities  ranging  from  1.1  X  10  to 
S  14  -3  . 

2.3  X  10  cm  This  data  is  used  in  Section  3.4.4  to  determine  the 
relevant  collision  cross  sections  controlling  the  effective  lifetimes  at 
the  higher  densities.  Before  proceeding  to  such  an  analysis,  however, 
we  must  first  establish  an  appropriate  method  for  calculating  the 
effective  radiative  lifetimes  of  the  two  upper  states  under  the  present 
conditions  where  the  resonance  radiation  from  each  state  only  escapes 
from  the  system  after  a  sequence  of  absorption-reemission  events.  This 
is  discussed  in  the  next  subsection. 


(a)  P8521  /  F8944 

Fit:  2.1% 


c  lit,- 8.9  ns 


T#=  31.6  ns 


40  60 

Time,  ns 


(b)  P8944/F8521 

Fit:  1.7% 


S  At,  =7.0  ns 


T  ,  -  30.5  ns 


Time,  ns 


Figure  3.12  (a)  t  (b) .  Exaaple  set  of  double  exponential  fits  to  the 

digitised  oscilloscope  traces,  taken  at  a  Cs  density  of 
1.36X10  cb  On  each  graph  the  points  represent  the 
digitised  data,  and  the  curve  is  the  theoretical  fit  having 
the  time  constants  indicated.  The  pump/fluorescence 
wavelengths  are  shown  in  each  case,  and  the  standard 
deviation  between  the  data  and  the  fit  is  given  as  a 
percentage . 
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Figure  3.12  (c)  k  (d) .  Example  set  of  double  exponential  fits  to  the 

digitized  oscilloscope  traces,  taken  at  a  Cs  density  of 
13  -3 

1.36X10  cm  .  On  each  graph  the  points  represent  the 
digitized  data,  and  the  curve  is  the  theoretical  fit  having 
the  time  constants  indicated.  The  pump/fluorescence 
wavelengths  are  shown  in  each  case,  and  the  standard 
deviation  between  the  data  and  the  fit  is  given  as  a 
percentage . 


3.4.3  The  Transport  of  Iesonance  Kadi at ion 

The  basic  theory  of  the  transport  of  imprisoned  resonance 
radiation  was  first  developed  independently  by  Holstein  (1947,  1951)  and 
by  Biberman  (1947) .  The  theory  has  subsequently  been  further  refined  by 
a  number  of  authors,  and  relevant  references  may  be  found  in  the  papers 
of  Irons  (1979  a,b,c).  Here  we  shall  draw  primarily  on  the  work  of 
Walsh  (1959). 

3.4.3. 1  The  Escape  Factor,  g 

It  was  shown  by  Holstein  (1947)  that  the  decay  of  the  density  of 
excited  states  within  an  enclosure  may  be  represented  by  the 
superposition  of  the  eigenmodes  characteristic  of  the  enclosure  shape. 

In  the  absence  of  an  external  source  of  further  excitation  the  spatial 
form  of  the  excited  state  density  rapidly  approaches  that  of  the  lowest 
eigenmode,  and  the  subsequent  local  rate  of  change  of  excited  state 
density  N*  may  be  represented  simply  by 


dn*/dt  =  -gN*/rn  (3.24) 

where  g  is  the  "escape  factor"  associated  with  the  lowest  order  decay 
mode.  The  value  of  g  depends  on  the  shape  of  the  container,  on  the  size 
of  the  container  measured  in  terms  of  the  optical  depth,  and  on  the  line 
shapes  for  emission  and  absorption.  These  quantities  depend  in  turn  on 
the  number  density  of  participating  atoms,  and  on  their  temperature. 
Asymptotic  expressions  for  g,  appropriate  for  large  optical  depths,  were 
calculated  approximately  by  Holstein  for  the  separate  cases  of  Doppler 
and  collision  broadening,  for  plane  parallel  slab  and  infinite  cylinder 
geometries.  Subsequent  more  detailed  calculations,  summarized  by  Irons 
(1979b)  have  shown  these  expressions  to  be  accurate  to  within  a  few 
percent.  In  the  present  problem  however  the  range  of  densities  covered 
is  such  that  we  must  in  general  take  into  account  the  combined  effects 
of  both  Doppler  and  collision  broadening,  and  in  addition  we  must 
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recognize  that  at  the  lowest  densities  the  asymptotic  form  for  g  will  be 
inaccurate,  and  must  in  some  way  be  modified. 

The  problem  of  including  the  effects  of  both  Doppler  and 
collision  broadening  in  the  calculation  of  g  was  addressed  by  Walsh 
(1959) .  Bis  analysis  led  him  to  propose  the  interpolation  formula 

g  =  gd  exp(-gcd2/gc2)  +  gc  E2(gcd/gc)  (3.25) 

where  Eg  is  the  error  integral.  The  subscripts  *c,"  "d,"  and  *cd*  refer 
respectively  to  the  escape  factors  applicable  to  a  collision  broadened 
line,  a  Doppler  broadened  line,  and  an  equivalent  escape  factor 
representing  the  effect  of  events  involving  Doppler  absorption  followed 
by  emission  in  the  relatively  transparent  collision  dominated  wings. 

Walsh’s  interpolation  formula  was  used  by  Zollweg,  Liebermann 
and  McLain  (1981)  in  a  slightly  modified  form.  These  authors  also 
introduced  modifications  to  the  asymptotic  forms  for  gd  and  gcd  such 
that  their  values  tended  to  unity  for  small  optical  depths.  Based  on 
this  work  we  have  used  the  following  formula: 


« =  gd  «P<-gcd2/gc2)  *  (gcd-3  *  scc~3r1/3 

(3.26) 

where 

gd  =  1.60/{l.«0  ♦  (k0B/fe)[t  ln(l  .  koB/fe)]1/2} 

(3.27) 

gc  -  «p(-k0R/feat1/2)  I0(k0»/fX/2> 

(3.28) 

gcd  '  */{a  *  (*/2)  [»  l»d  *  k0B/fe)J1/2} 

(3.29) 

a  =  (l/Nrc)(XoN/4r)(M/2kTc)1/2 


(3.30) 


ko  *  (X03N/8i3/2) (,u/tl) (l/ro) (H/2kTe) 


(3.31) 


These  equations  differ  from  those  of  Zollweg  et  al.  (1081)  in 
three  respects.  The  najor  change  is  that  the  effects  of  natural 
broadening  are  oaitted  froa  the  definition  of  the  Voigt  parameter  *&", 
proportional  to  the  ratio  of  the  collision  broadened  width  to  the 
Doppler  width.  This  change  is  largely  pragmatic,  in  that  we  then  find 
the  interpolation  foraula,  Eq.  (3.26)  or  Eq.  (3.25)  to  be  better  behaved 
in  the  region  kQR  -  1  and  lower,  and  in  the  transition  between  the 
regions  of  predominantly  Doppler  broadening  and  collisional  broadening. 
The  omission  of  natural  broadening  effects  can  be  rationalised  on  the 
basis  of  the  anticipated  incomplete  frequency  redistribution  in  the 
natural  wing,  as  discussed  by  Huennekens  and  Gallagher  (1083) . 

The  interpolation  formula,  Eq.  (3.25)  or  Eq.  (3.26),  is  the 
original  fora  proposed  by  Walsh  (1959) .  We  find  the  modification  made 
by  Zollweg  et  al .  (1981) ,  to  control  the  asymptotic  value  of  g  as  kQR 
becomes  very  small,  to  be  unnecessary  provided  we  use  the  modified 
definition  of  the  Voigt  parameter  given  in  Eq.  (3.30). 

The  remaining  difference  lies  in  the  constant  chosen  for 

Eq.  (3.27).  Zollweg  et  al.  (1981)  used  a  value  of  1.82  for  consistency 

with  the  results  of  Phelps  and  McCoubrey  (1960)  for  5  <  kQR  <  400,  which 

are  approximately  20X  larger  than  are  given  by  Holstein’s  asymptotic 

foraula.  We  find,  however,  that  use  of  the  value  1.60  gives  better 

agreement  at  low  kQR  values  with  the  results  of  Blickensderfer  et  al. 

(1976) ,  and  clearly  will  give  the  correct  asymptotic  limit  at  large 

values  of  k  R. 

o 

Evaluation  of  Eq.  (3.26),  and  the  associated  Eqs.  (3.27)  through 
(3.31),  proceeds  separately  for  the  two  resonance  lines  involved  in  the 
present  context.  For  the  natural  lifetimes  we  use  the  values 
rn(6%*^2)  =  34  ns  and  =  30-5  n8  1988;  Gallagher, 

1967) .  For  the  density  normalised  collisional  broadening  frequencies  we 
use  the  values  (l/Nf^jyg  =  8.8  x  10~®  s  *  and  (l/NT^gyg  = 
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1.51  x  10  s  ,  given  by  Chen  and  Phelps  (1968).  The  values  used  for 
f  ,  the  effective  number  of  hyper fine  components,  were  deduced  from 
considerations  of  the  ratio  of  the  separations  of  the  components  to  the 
trapped  width  of  the  individual  components,  as  discussed  in  the  next 
subsection. 

3. 4. 3. 2  The  Effective  Number  of  Hyperfine  Components,  f 

In  many  instances  the  radiation  of  a  given  line  consists  of  a 
number  of  closely  spaced  but  distinct  components  due  to  the  presence  of 
various  isotopes  and/or  hyperfine  splitting  of  the  level  arising  from 
the  interaction  of  the  orbital  electrons  with  the  nuclear  spin.  Cesium 
has  only  one  stable  isotope  so  we  are  concerned  here  only  with  the 
latter  effect. 

In  the  simplest  possible  case,  where  the  line  consists  of  "f* 
equally  intense  components,  separated  sufficiently  from  each  other  that 
the  levels  do  not  exchange  energy  radiatively,  the  effect  of  the 
hyperfine  structure  on  the  radiative  transport  is  accounted  for  by 
simply  dividing  the  optical  depth  kQR  by  the  number  of  components,  f. 

In  effect,  each  line  component  "sees*  an  absorber  density  of  N/f.  At 
the  other  extreme,  where  the  separation  of  the  components  is  much 
smaller  than  the  line  widths,  the  line  structure  has  no  effect  on  the 
process  of  radiation  transport.  In  general,  however,  the  situation  is 
more  complicated,  due  to  the  partial  overlap  of  some  or  all  of  the 
components,  and  also  to  a  lesser  degree  by  the  fact  that  different 
components  have  different  relative  intensities.  For  present  purposes  we 
take  account  of  this  complication  through  the  use  of  an  effective  number 
of  components,  represented  by  the  factor  f  in  the  equations  given  in 
the  previous  subsection.  To  estimate  f  we  use  the  same  general 
approach  as  was  employed  by  Walsh  (1959) ,  without  the  level  of  detailed 
calculation  which  he  performed  for  mercury. 

Following  Holstein  (1947) ,  we  note  that  for  a  multicomponent 
line  the  net  transmission  coefficient  may  be  written 


mm  BWl^roWTOfTOOffiCTWB 
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T(/>)  =  E  x.  T.(/>)  (3.32) 

where  p  is  distance,  is  the  intensity  of  the  i’th  component,  and  E 
x^  =  1.  I\  (/>)  is  the  transmission  coefficient  for  the  i’th  component, 
for  which 


k. 

1 


=  x. 
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k 

o 


(3.33) 


For  the  purpose  of  this  discussion  we  assume  sufficiently  large 
values  of  k ^p  at  line  center  for  the  asymptotic  forms  given  by  Holstein 
to  apply  for  each  isolated  component.  For  example,  when  Doppler 
broadening  dominates 


T.(/0  =  (l/k./0(l/flnk./01/2  (3-34) 

To  estimate  T^(p)  for  a  particular  component  in  the  presence  of  adjacent 
components  we  must  consider  the  extent  to  which  the  lines  overlap.  The 
escape  of  radiation  from  a  particular  line  occurs  at  wavelengths 
sufficiently  displaced  from  line  center  that  k^  i  1.  Other 
wavelengths,  within  the  absorbing  core  are  essentially  trapped. 

Following  Valsh  (1959)  we  define  an  absorption  half-width  for  each 
component  by 


kfO'o  *  *i>  =  1  (3 •  33) 

For  a  Doppler  broadened  line  we  obtain 

*d  =  (Vxo><lnV/f>1/2  (3-36) 

while  for  a  collision  broadened  line  we  have 
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Ve  nay  discuss  the  degree  of  overlap  of  lines  separated  by  A  in 
terms  of  the  ratio  A/f. 

Lines  Completely  Isolated.  L21  * 

Under  these  conditions  use  of  Eq.  (3.32)  is  straightforward, 
giving  in  the  case  of  Doppler  broadening 

Td(/0  =  E  [l//>ko(Tlnx.ko/»)1/2]  (3.38) 

This  may  be  written  in  the  form 


T dW  =  C(,lko/f)'1[rl„(/,kt)/()]-1/2  (3.39) 

•her.  C  =  r1[lo(^ko/f)]1''2  E[ln(/>xik0)]1/2  (3.40) 

is  close  to  unity,  and  is  identically  so  when  the  f  components  are  of 
equal  intensity,  i.e.  when  x^  =  f  all  i. 

Similarly,  for  isolated  collision  broadened  components  we  obtain 

T c(/0  =  [a/r1/2/,(ko/f)]1/2  IKx./f)1/2  (3.41) 

where  the  sum  is  close  to  unity  and  again  is  identically  so  when  the  f 
components  are  of  equal  intensity.  Hence,  to  a  good  approximation,  the 
transmission  coefficient  for  a  line  consisting  of  f  isolated  components 
is  obtained  from  the  appropriate  form  for  a  single  line  by  using  kQ/f  in 
place  of  kQ.  Therefore,  in  this  situation  we  assume  the  effective 
number  of  components,  f  ,  to  be  equal  to  the  actual  number  of 
components ,  f . 


T(/>)  =  \  [xjTjC/j)  +  xfTfOO]  (3.42) 

where  T^(p)  and  are  the  values  which  would  apply  to  these 

components  if  they  were  isolated.  For  Doppler  broadened  lines  these  are 
given  by  Eq.  (3.34),  and  hence  we  obtain  from  Eq.  (3.42) 


Tdw  =  01/[^ko(.ln,.ko)1/2] 

(3.43) 

where 

ci  =  2  „>1/2  (a°vk0>‘1/2  *  u»yk0r1/2] 

(3.44) 

=  [lnA  /ln(,k  /f)]l/2 

Thus  the  factor  is  in  general  >  1,  and  tends  to  unity  for 

large  />kQ,  when  Eq.  (3.43)  gives  the  value  appropriate  for  complete 
overlap  of  all  the  components.  Therefore  in  this  situation  we  expect 
the  effective  number  of  components,  f  ,  to  be  in  the  range  1  if  <  f. 

For  collision  broadened  lines  Eq.  (3.42)  gives 


T.00  =  [a/T1/2,>k0] 1/2  [|  (*!1/2  ♦  x//2)] 


(3.45) 


=  [a/ir*/2pkof]  */2  when  x^  =  f  *  ,  all  i. 


We  note  that,  in  this  region  (i.e.  partial  overlap  of  collision 
broadened  lines)  the  transmission  coefficient  is  lower  than  the  value  it 
would  attain  with  complete  overlap,  approximately  by  the  factor  f  ^2. 


The  effective  number  of  components,  f  ,  is  therefore  predicted  to  go 

_1  ® 

through  a  minimum  value  of  f  before  approaching  unity  when  the  lines 
completely  overlap. 

The  hyperfine  structure  relevant  to  the  present  problem  is  shown 

♦  ♦ 

in  Figure  3.13  for  the  8944  A  line,  and  in  Figure  3.14  for  the  8521  A 
line.  The  hyperfine  energies  have  been  calculated  from  the  standard 
equation  (Arimondo  et  al.  1977,  Eq.  2.7) 

T(F)  =  \  AK  +  B[(3/2)K(K+1)  -  21(1+1) J(J+1)]/2I(2I-1)2J(J-1)  (3.46) 

where  K  =  F(F+1)  -  1(1+1)  -  J(J+1) 

using  the  values  of  A  and  B  recommended  for  Cs  in  the  same  reference, 
namely 

A(6pi/2)  =  291.90  MHz;  A(6p3/2)  =  50.34  MHz;  B(6p3/2)  =  -0.38  MHz 

In  Figures  3.13  and  3.14  the  fractional  relative  intensities  are 
also  indicated.  The  values  shown  were  obtained  from  the  tables  of  White 
and  Eliason  (1933) . 

The  values  of  ♦  ,  and  #  have  been  calculated  using  Eqs.  (3.36) 

u  C 

and  (3.37)  for  each  set  of  experimental  conditions.  The  values  of  •j, 

calculated  using  the  actual  values  of  f,  range  from  251  to  783  MHz  for 

the  8944  A  line,  and  from  291  to  841  MHz  for  the  8521  A  line.  We  note 

♦ 

that  the  components  of  the  8944  A  line  begin  to  overlap  due  to  Doppler 
broadening  only  at  the  highest  densities.  Therefore,  in  computing  gj 
from  Eq.  (3.27),  and  gcd  -from  Eq.  (3.29) 7  we  have  assum«d  'a  Value  of  -  -  - 
fe(8944)  =  4  at  low  densities,  and  reduced  its  value  at  higher  densities 
according  to  the  algorithm 

f e(8944)  =4/(l  +  *d (8944) /I 168) 


(3.47) 
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Figure  3.13  The  hyperfine  levels  ^>f_  the  ^  end  states  of  Cs 

and  the  resulting  hyperfine  structure  cf  the  8944  A  line. 


The  separation  and  relative  weights  of  the  four  components 
are  shown  in  the  lower  right  portion  of  the  figure. 
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In  computing  •  (8944)  froa  Eq.  (3.37)  we  have  used  f  =4 
throughout.  The  resulting  values  range  from  0.6  to  1201  MHz.  When 
evaluating  g  (8944)  from  Eq.  (3.28),  we  have  used  for  f  (8944)  the  value 

C  6 

4  when  ♦  (8944)  <  584  MHz,  and  the  value  2  otherwise. 

c  • 

The  six  components  of  the  8521  A  line  fora  two  distinct  groups, 

and  the  three  components  of  each  group  overlap  even  at  the  lowest 

density.  Therefore,  in  computing  g  (8521)  and  V  (8521)  we  have  used  a 

c  c 

value  of  f  (8521)  =  2  throughout.  The  values  of  •  (8521)  computed  from 
c  c 

Eq.  (3.37)  range  froa  1.6  to  3231  MHz.  Hence,  as  the  density  is 

increased  fe(8521)  will  range  froa  6  to  approximately  1/2,  the  latter 

number  arising  from  the  partial  overlap  of  two  lines,  each  of  which 

consists  of  three  completely  overlapping  components.  However,  over  the 

range  where  collision  broadening  is  important  the  effective  value  of 

fe(8521)  is  2,  and  for  simplicity  we  have  assumed  this  value  throughout. 

These  choices  of  f  are  summarized  in  Table  3.3. 

© 

Table  3.3  Showing  the  values  chosen  for  f  ,  the  effective 
number  of  hyperfine  components,  in  calculating 
the  quantities  indicated.  See  Figures  3.13  and 
3.14,  and  text  for  discussion. 


Values 

of  f 

used  to 

calculate : 

Line 

f 

*d 

* 

c 

*d 

gcd 

gc 

8521 

6 

6 

2 

2 

2 

2 

8944 

4 

4 

4 

4*2 

4*2 

4*2 

f 


8.4.4  Keconciling  Computed  and  Measured  Values  of  and 

Values  of  and  Tg  are  calculated  using  Eqs.  (3.21)  and  (3.22) 
respectively.  In  so  doing  the  value  of  is  computed  from  the  equation 


K1  -  'l/Tnl  *  'T^ql 


(3.«) 
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where  is  coaputed  as  discussed  in  the  preceding  subsection, 
the  self  quenching  cross  section  for  level  1,  treated  here  as  an 
adjustable  input  parameter,  and  v,p  is  the  average  center-of-aass  thermal 
velocity,  given  by 


vT  =  4(kT./wM)1/2  (3.49) 

A  similar  equation  to  Eq.  (3.48)  applies  for  computation  of  K^,  with  a  g 
being  treated  as  an  adjustable  input  parameter.  The  values  of  ud 
K 2i  are  coaputed  by  assuming  a  trial  value  for  the  exothermic  cross 
section,  and  calculating  the  endothermic  cross  section  from  detailed 
balancing.  For  example,  in  the  case  where  u ^  we  assume  a 

temperature  independent  trial  value  for  a ^  end  use 

a2\f°\2  =  ^ (36°) 

to  calculate  the  corresponding  trial  value  of  Multiplying  these 

cross  sections  by  Nv^  then  gives  and  respectively.  The  Cesium 
atomic  density  N  is  given  by  the  relationship  (Taylor  and  Langmuir, 

1937) 


log10[ffTc/9.668  *  1018]  =  11.0631  -(4041/Tr)  -  log1()Tr  (3.51) 

In  order  to  j  udge  the  overall  quality  of  the  fit  of  the  theory 
to  the  experiment,  and  as  an  aid  in  optimising  the  choice  of  the  fitting 
parameters  a ^ ^ ■  and  a we  used  commercially  available  software  to 
compute  and  screen-plot  the  values  of  and  as  functions  of  the 
computed  Cesium  density,  together  with  the  measured  values  determined  by 
the  method  described  in  subsection  3.4.2  This  allowed  trial  and-error 
fitting  to  proceed  rapidly. 


It  was  found  possible  to  fit  the  measured  values  of  r.  and  r 

-13  2  is 

with  any  value  of  a19  £  1.4  x  10  cm  ,  by  appropriately  adjusting  a  . 

* ®  _i3qi 

and  a  The  limiting  sets  of  cross  sections  are:  (a10  =  1.4  x  10 
2  v*  -14  2  _  i  *3 

cmz ,  .  =  0,  and  a  9  =  3  x  10  cm  j ,  and  (ff10  =  0,  a  ,  =  2  x  10 

2  V  ”  -14  2  “  V 

cm  ,  and  a  2  =  2  x  10  cm  ) .  In  order  to  further  constrain  the  range 

of  acceptable  values  we  have  also  attempted  to  fit  the  measured  relative 

signal  strengths  of  the  direct  and  transferred  fluorescence  signals. 

For  convenience  we  choose  to  do  this  at  the  instant  of  the  peak  in  the 

transferred  fluorescence  signal,  given  by  Rq.  3.23.  The  ratio  of  the 

signal  strengths  is  given  by 


Where 


^2  _  X1  Tnl  *2 


»i  ( v 


n2 


(3.52) 


Ni(tp)  -  K12 


««P  C-tp/rB) -exp (~tp/rf) 


X)  exp  (~t  /rf )  -  (I1+K12-rf _1)  exp (-t  JrJ 


(3.53) 


is  the  ratio  of  the  densities  of  the  two  excited  states,  and  g,  0  are 
calculated  as  discussed  in  Section  3.4.3.  These  equations  apply  to  the 
case  where  level  1  is  pumped  and  the  signal  from  level  2  peaks  at  time 
tp.  The  corresponding  equations  for  the  case  where  level  2  is  pumped 
are  obtained  by  simply  interchanging  the  subscripts.  The  signal  ratios 

predicted  in  this  way  are  compared  to  the  measured  ratios  corrected  for 

the  wavelength  dependence  of  the  photomultiplier  sensitivity,  and  the 
differing  transmissions  of  the  interference  filters  used  to  isolate  the 
recorded  signals. 

The  best  fit  to  the  measured  signal  ratios  is  obtained  with 
ffj 2  =  1-3  x  10  ^  ca^,  and  =  ff^2  =  2  x  10  ^  crn^,  shown  in  Fig. 
3.15.  The  fit  to  the  measured  values  of  and  Tg  achieved  with  these 
values  of  the  cross  sections  is  shown  in  Fig.  3.16.  This  cross  section 

set  gives  the  best  overall  fit  to  the  present  data. 
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Figure  3.15  The  ratio  of  the  fluorescence  signals  at  the  instant  t^  of 
the  peak  of  the  transferred  signal.  Subscripts  1  and  2 
refer  to  the  upper  (($^-3/2)  and  the  lower  (S2!^^)  excited 
states  respectively.  Points  are  experimental  data  and  the 
full  curves  are  predictions  of  the  model  using  '12  =  13  X 
10~13  cm2  and  =  a q2  =  2  x  10  1  cm2. 


A  previous  measurement  of  a. _  by  Cjakowski  and  Krause  (1965), 


performed  at  auch  lower  densities  where  the  effects  of  radiation 


15  o 

trapping  could  be  ignored,  gave  a  value  of  3.1  x  10  cm  .  This  value, 


—13  2  —14  2 

together  with  <7^  =  2  x  10  ca  and  =  2  x  10  cm  ,  provides  an 


acceptable  fit  to  the  present  and  data,  shown  in  Fig.  3.17. 
However,  this  set  of  cross  sections  fails  significantly  in  fitting  the 


signal  ratio  data,  as  shown  in  Fig.  3.18.  Thus,  we  are  forced  to  the 
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conclusion  that  the  value  of  =  3.1  x  10  ca*  is  incoapatible  with 


the  present  measurements .  No  explanation  for  this  discrepancy  with  the 
results  of  Cjakowski  and  Krause  (1965)  can  be  offered  at  this  time. 

For  comparison  purposes  we  also  show  in  Figures  3.18  and  3.17 
the  behavior  predicted  for  the  case  of  sero  quenching  cross  sections. 

The  maximum  value  of  r&  achieved  in  this  situation,  in  the  present  size 
of  container  (R  =  5.4  ca) ,  would  be  83  /is,  compared  to  35  /is  achieved  in 
the  presence  of  the  rather  large  quenching  cross  sections  needed  to 
account  for  the  data. 


3.5  PUMP  P0IB1  IMPLICATIONS 


The  conclusions  drawn  above  regarding  the  achievable  effective 
lifetime  of  the  Cs*(8^P)  states  allow  us  to  Bake  a  firmer  estimate  than 
was  previously  possible  of  the  puap  power  density  likely  to  be  needed  in 
an  optimized  Cs-Tl  ALF  module.  For  example,  assuming  the  module  has  the 
same  size  as  the  experimental  tube  used  for  the  Cs  fluorescence  work 


reported  above,  and,  as  in  the  calculations  of  Section  2,  that  we  need 


10  -3 

to  sustain  a  T1  density  of  10  ca  ,  and  that  this  would  involve  a 


10  -3 

corresponding  density  of  4  x  10  ca  of  T1q,  we  can  use  Bq.  (3.0)  to 


calculate  the  dependence  of  [Cs*]  on  [CsQ] .  The  implied  resonance  line 
puap  power  density  required  to  sustain  the  Cs*  is  then  calculated  using 
appropriate  to  the  value  of  [Csq] ,  read  from  rg  of  Fig.  3.15.  The 
resulting  required  pump  power  density  has  a  minimum  value  of 


3  13  -3 

approximately  6  ml /cm  in  the  region  where  [C»  ]  =  10  ca  ,  with  a 

^  11  o  O 


*  11-3 

corresponding  value  of  [Cs  ]  =8  X  10  cm  These  densities  of  Cs  and 
T1  are  achievable  simultaneously  at  a  temperature  close  to  400*C 
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3. lx  10" 40  cm 


Figure  3.18  3 mm  as  Fig.  3.15  except  that  a.„ 

-13  2  **  -14  2 

8  .  «  2  x  10  cm  ,  Mid  8  a  ■  3  x  10  ca  have  been  used 
ql  q* 

to  coapute  the  full  curves. 


ra 


i  .  1 


provided  one  uses  an  amalgam  of  the  appropriate  composition  (Hirayama, 
1986) . 

The  above  estimate  applies  to  the  locally  absorbed  pump  power 
density,  and  therefore  represents  a  lower  bound  on  the  resonance  line 
power  that  would  need  to  be  applied  to  a  given  module.  Because  the 
medium  is  optically  thick  it  will  be  desirable  to  employ  a  Cs  resonance 
line  source  giving  a  suitably  self-reversed  profile,  such  that  the 
module  has  an  effective  optical  thickness  of  approximately  unity.  In 
this  way  it  should  be  possible  to  efficiently  pump  the  module  more  or 
less  uniformly.  Provided  the  pump  line  profile  can  be  tailored  in  this 
way  the  actual  supplied  optical  power  would  probably  only  need  to  be  a 
factor  of  two  or  three  greater  than  the  lower  bound  obtained  from  the 
above  calculation. 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  section  we  present  separately  the  conclusions  and 
recommendations  which  follow  from  the  work  described  in 
Section  2,  concerning  the  T1C1  based  ALP  concept,  and  from 
the  work  described  in  Section  3  which  is  pertinent  to  the 
Cs-Tl  based  ALP. 

4.1  THB  T1C1  BASED  ALP  CONCEPT 

The  following  conclusions  are  drawn  from  the  sample 
calculations  presented  in  Section  2: 

1.  There  exists  an  extensive  region  of  parameter  space  where  a 

10  -3 

density  of  approximately  1  X  10  cm  of  thallium 

metastables  can  be  maintained  by  photodissociation  of  T1C1 

2 

using  pump  power  of  approximately  100  mV/cm  .  The  latter 
number  is  approximately  inversely  proportional  to  ,  the 
cross  section  for  production  of  T1  by  photodissociation  of 
T1C1.  " 

2.  Mercury  resonance  line  radiation  at  254  nm  is  suitable  for 
pumping  this  particular  process. 

3.  The  required  densities  of  T1C1  are  sufficiently  low, 
typically  1  x  10*^  ca  that  the  system  is  optically  thin 
to  the  pump  radiation.  This  allows  large  volumes  to  be 
pumped  uniformly  without  resorting  to  specially  tailored 
pump  line  profiles  or  geometries. 

4.  The  required  container  temperature  has  a  lower  bound  of 
approximately  250*C,  set  by  the  prescribed  density  of  T1C1, 
but  a  somewhat  higher  temperature  may  be  necessary  to  avoid 
excessive  coverage  of  the  windows  by  condensed  thallium. 


The  actual  value  could  best  be  determined  by  experiment,  but 
is  not  expected  to  exceed  400* C. 

5.  Tbe  density  of  T1q  is  conveniently  controlled  by  the 

addition  of  excess  chlorine  to  tbe  mixture.  An  amount 

15  “3 

sufficient  to  give  [Cl]  SIX  10  cm  ,  corresponding  to 
12  -3 

[Clg]  -  1  X10  cm  ,  is  suitable  for  this  purpose. 

6.  With  the  values  assumed  for  the  buffer  gas  quenching  cross 

-19  2 

section  (oa  =  1  X  10  cm  )  and  the  three  body 

o  —32  6 

recombination  rate  coefficients  (kg  =  =  1  X  10  cm  /s) 

a  buffer  gas  cannot  be  added  in  sufficient  amount 
18  -3 

(>1  X  10  cm  )  to  usefully  promote  the  volume 

recombination  processes.  Smaller  amounts 
16  —3 

([M]  -  3  x  10  cm  )  enhance  the  sustained  density  of 
metastables  by  inhibiting  their  transport  to  the  malls,  but 
at  the  expense  of  also  increasing  [T1q] . 

7.  Because  of  the  direct  sensitivity  of  the  required  pump  power 
to  the  value  of  (the  cross  section  for  Tla  production  by 
photodissociation  of  T1C1  at  254  ns)  and  in  view  of  the 
uncertainties  present  in  the  literature  regarding  its  value, 
an  independent  measurement  of  this  quantity,  and  if  possible 
of  on  (the  corresponding  cross  section  for  T1  production) , 
would  be  very  valuable. 

8.  This  concept  continues  to  appear  viable,  but 

proof -of -principle  experiments  using  CW  pumping  are  needed 
to  better  establish  the  required  operating  temperature  and 
pump  power. 

4.2  THB  Cs-Tl  BASBD  ALT  CONCEPT 

The  following  conclusions  and  recommendations  are  based  on  the 
work  discussed  in  Section  3  of  this  report. 

1.  The  measured  effective  lifetimes  and  relative  fluorescence 
signals  of  the  Cs(6% jyg)  *od  Cs(6^P  <g)  in  pure  Cs  vapor 
are  consistent  with  lowest  order  mode  radiation  trapping 


modified  by  collisional  exchange  of  radiation  between  the 

-13  2 

states  with  an  exothermic  cross  section  of  1.3  X  10  cm  , 

-14  2 

and  a  self  quenching  cross  section  of  2  X  10  cm  of  both 
states  by  ground  state  Cs. 

The  effect  of  the  collisional  self  quenching  is  to  define  an 
optimum  Cs  density  at  which  the  effective  lifetime 
maximises,  at  ~35  fia  in  the  present  apparatus.  In  the 
absence  of  quenching  a  lifetime  of  83  /la  could  have  been 
achieved.  The  lower  value  implies  higher  pump  power  than 
would  otherwise  be  needed,  but  is  nevertheless  acceptable  in 
the  present  context. 

For  a  module  of  the  same  dimensions  as  the  present 

experimental  tube,  operating  at  approximately  400* C,  the 

required  applied  optical  pump  power  would  be  approximately 
3 

10-20  mV/cm  provided  the  pumping  line  profile  and  the  Cs-Tl 
amalgam  composition  were  both  optimised. 

A  possible  source  of  background  noise  in  a  Cs-Tl  ALF  is 
in-band  radiation  from  higher  lying  Cs  states,  such  as 
7^Pjy2  g/g  *  giving  radiation  at  450  and  455  nm 

respectively.  The  mechanism  of  excitation  of  such  states 
would  of  necessity  be  a  second  order  effect  (Leslie  et  al., 
1977) ,  but  the  resulting  signal  might  nevertheless  be 
significant.  The  Cs  cell  used  for  the  work  described  in 
Section  3  provides  the  opportunity  to  observe  such  radiation 
in  pure  Cs  following  resonance  line  pumping.  Measurements 
of  this  type  have  been  performed  under  another  program,  and 
will  be  reported  in  detail  elsewhere  (Liu  et  al.,  1086). 
Based  on  the  measured  signal  dependences  on  the  Cs  density, 
and  on  the  pumping  power  density,  the  observed  signals  are 
ascribed  to  two-photon  absorption  of  CsQ  at  the  high 
instantaneous  pump  power  densities  used  in  the  measurements. 
Under  the  low  power  Cf  pumping  conditions  prescribed  for  an 
actual  ALF  the  resulting  background  is  predicted  to  be 
negligibly  small. 


As  the  next  major  step  in  demonstrating  the  viability  of  a 
Cs-Tl  ALF  work  is  in  progress  using  a  Tl-Cs  cell.  The  first 
objective,  to  observe  the  emission  of  378  nm  radiation 
resulting  from  the  sequence  of  events  represented  in 
Figure  3.1,  has  recently  been  achieved.  A  detailed  account 
of  this  work  will  be  reported  in  detail  elsewhere 
(Liu  et  al. ,  1986) . 

The  remaining  uncertainties  in  determining  the  viability  of 
this  concept  are  the  cross  sections  appearing  in  Eq.  (3.0). 
Of  most  importance  is  the  ratio  of  a  determination  of 

which  can  best  be  made  using  a  Cs-Tl  cell  in  which  the 
concentrations  of  T1q,  Tl^,  Csq,  and  Cs*  can  be  individually 
monitored  as  a  function  of  the  optical  pumping  of  the  Cs. 

Ve  plan  to  perform  measurements  of  this  type  in  the  near 
future.  Knowledge  of  all  three  cross  sections,  5  6’  *or 
reactions  2,  5,  and  6  in  Table  3.1,  is  a  prerequisite  if 
confidence  is  to  be  placed  in  the  design  of  a  prototype  ALF 
module. 

Based  on  our  best  estimates  of  these  unmeasured  cross 
sections,  we  are  unaware  of  any  factor  that  would  prevent 
the  demonstration  of  a  practical  Tl-Cs  ALF. 
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